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Abstract: 


Dystric/eutric Cambisols diagnostic features from TEFR Petrohan were investigated 
in the present study, which morphologically corresponded to Podzols. Diagnostic horizons 
were investigated to establish presence or absence of podsolization process and to classify the 
soil type according to the modern requirements of WRB, (2006, 2007; 2014). Humus 
forming process insitu, clearly expressed in the A, horizon, acidification, degradation of clay 
minerals, formation of EF horizon, which is clay-poor with increased content of sand, 
accumulation of Fe and other mobile elements in Bw horizon were established. The 
metamorphic Bw horizon had some diagnostic properties, which met these of spodic horizon. 
This horizon is specific for the Podzols. The content of organic carbon in it was less than 0.5 
%. Cambisols properties are predominant.It was established that the soil-forming process in 
the investigated Cambisols is directed to podsolization process. 


Key words: Cambisols, Podzols, org.C, albic, spodic diagnostic horizon 


Introduction 

Dystric/eutric Cambisols are significant for the forestry in Bulgaria. Depending on the 
combination of the main soil forming factors, in some of them are established specific 
processes. For example, one of them is the podsolization process that creates conditions for 
the formation of a new soil type called Podzols. 

Some of the first sources of information for the Podzols are mentioned in Russian 
literature in the 1950s (Williams, 1949, Garkusha, 1951, Vilenskii, 1954, Russel, 1955, 
Alexandrova et al., 1958). In the following years, a number of authors (Petersen, 1976, 
Mokama and Buurman, 1980; Anderson et al., 1982; Petersen, 1984; Buurman, 1984; 
Lundstrém et al., 2000) studied their genesis, chemical, mineral composition and key 
features. 

According to the World Reference Base for Soil Resources (WRB, 2006, 2007) and 
the International Union of Soil Sciences (USS) Working Group WRB, 2014 Podzols covers 
485 million ha worldwide, mainly in the temperate and boreal zone of the northern 
hemisphere. They occupy vast territories in Scandinavia, Northwest Russia and Canada. 

The first information in Bulgaria about these soils is given by Boyadzhiev in 1967. It 
is known that they are distributed mainly in the middle and high mountain areas. They are 
formed under the influence of coniferous and mixed forests, humid, cold climates and silicate 
rocks. According to Achkov's studies (1980, 1981), the soils that correspond to in the country 
are Cambisols. They can evolve at a certain stage of their development into Podzols. This 


Bulgarian Journal of Soil Science® 2018 Volume 3. Issue 2. www.bsss.bg sl 


Petrova, 2018 





statement is confirmed by other authors (Dialo, 1995; Teoharov et al., 1995). Cambisols can 
change and form a spodic diagnostic horizon, which is specific for the Podzols. In them it is 
formed under the influence of climate change and the dominant vegetation, anthropogenic 
impact, etc., which alter the processes of decomposition of the litter and lead to the 
accelerated formation of water-soluble organic acids (Lundstrém et al., 2000; Sauer et al., 
2007). 

In Bulgaria a podzolization process of Umbrisols was described by Achkov (1980) in 
the high parts of the Rila mountains - 2485 m and an initial stage of podzolization process at 
an altitude of 1460 m - 1470 m also in Rila (1981). The presence of a typical spodic 
diagnostic horizon is mentioned by Dialo (1995), at an altitude of 1450 m in Pirin mountains. 
According to Dialo (1995), these soils are still not well described in the country, but they are 
a separate taxonomic unit and should find a place in the soil classification of Bulgaria. 
Boyadjiev has also mentioned the Podzols (1994a, 1994b), when he has established a 
correlation between the American Taxonomic System and that of the revised FAO- 
UNESCO-ISRIC legend with the Basic Soil Classification in Bulgaria. It is clear that the 
correlation is completely possible and necessary. 

According to unpublished data of Yordanov (1985) on the territory of the TEFR 
Petrohan there is a Podzols with clear morphological features. For this soil, no parameters 
have been studied so far that characterize the elementary soil-forming processes and 
diagnostic features to prove the albic and spodic horizons necessary for the definition of a 
Podzols, according to WRB requirements (2006, 2007 and 2014). The aim of the present 
study is to investigate the diagnostic features of Cambisols from the territory of the TEFR 
Petrohan, which morphologically corresponds to Podzols, to determine the diagnostic 
horizons, to establish the presence or absence of podzolization process and to classify the 
soil type according to WRB's modern requirements (2006, 2007, 2014). 


Materials and Methods 

The investigated soil is located at 1400 m altitude at the southern exposure. It has a 
depth of 30 cm and more. The soil-forming rock is granite. The litter is 1 cm, torn by grass 
and shrubs mainly covered by blackberry (Vaccinium myrtillus L.) and blueberry (Vaccinium 
uliginosum L.). When plotting the soil profile, morphological features different of these of 
Cambisols and typical for Podzols were established. The strongest impression was the colour 
change of separate horizons that point to the presence of albic and spodic diagnostic horizons. 
In the soil profile are defined AturfAEBC horizons by morphological features. 

Parameters that characterize diagnostic albic and spodic horizons are investigated: 
colour by horizons (Munsell scale); Texture - Modified Kachinski method (dispersion with 
sodium pyrophosphate); pH (CaCl12) - ISO 10390; pH (H20) - ISO 10390; Organic carbon - 
Modified Turin’s method (Kononova, 1963; Filcheva E., S. Tsadilas, 2002); exchange cations 
- ISO 11260 & ISO 14254, 0.1 M / 1 BaCl2 solution and AAS Perkin Elmer 5000; 
exchangeable acidity - ISO 11260 & ISO 14254 (extract with 0.1 mol / 1 BaCl» solution to 
equilibrium desorption and titration with 0.05 mol / 1 NaOH); Fe content - by decomposition 
with Aqua regia and Flame AAS (ISO 11466) Perkin Elmer 5000; cation exchange capacity - 
determined as a sum of the basic cations and the exchangeable acidity. 
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Results and Discussion 

Studies on the E albic diagnostic horizon showed full compliance with WRB 
requirements (2006, 2007 and 2014). It laid directly under the humus-accumulative horizon. 
The distinctness between them was clear. The topography between the E and the B horizon 
was wavy with slightly albeluvic features. The soil’s 
colour in the EF horizon was bright - 10 YR 6/1 (when 
dry). It met the requirement for a value of less than 7 


(Table 1). 


Table 1. Soil colour by horizons (Munsell scale) 


Horizon | Soil Distinctness | Colour 
Depth 
(cm) 


CC 
[as__[07 [ear [YRS 





Photo 1. Soil Profile 


On the field, the Bw horizon was determined as a spodic diagnostic horizon. The 
requirements of WRB (2006, 2007; 2014) for its definition are: for colour, soil depth, amount 
of org. C, pHy20 and absence of a natric diagnostic horizon, with a high content of 
exchangeable Na, and / or Mg. 

The study showed that the soil colour in the Bw horizon corresponded to the 
requirements. It was determined as 7.5 YR 4/4 (when moist) and met the requirement for hue 
7.5YR, value 5 or less and chroma 4 or less. It laid directly below the albic diagnostic 
horizon. The soil depth of Bw horizon was 10 cm, which met the requirements for spodic 
horizon - minimum 2.5 cm. The pH of the soil solution was - pHy20 = 5.2 and also met the 
requirements for a spodic diagnostic horizon - pHy 0 of less than 5.8 (Table 2). 

A disparity with the spodic diagnostic horizon requirements of the WRB (2006, 2007; 
2014) was established only with respect to the amount of org. C, which was very low and did 
not exceed the minimum of 0.5% . The result was obtained as an average sample of the entire 
10 cm soil depth of the horizon (Figure 1). 
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Content of org. C, (%) 


A o7 


A 7-14 


E 1420 0.72, 


Bw 20-30 0/32 


Cc 30-\0li4 
0 2 4 6 8 10 
Figure 1. Content of org. C, (%) in the soil profile 


In order to trace the distribution of org. C in Bw horizon, it was divided into 4 parts by 
2.5 — as the required thickness for spodic diagnostic horizon (Figure 2). 


Content of org. C, % 


AN 0 Se )=4 12 
Ab? 5.10 Ds 2.24 
A3 10-146 ee 1s 
E1420 MEE 0.76 
Bw 20-22.5 MM 0.26 
Bw2 22.525 MMM 0.31 
Bw3 25-27.5 ME 0.30 
Bw4 27.530 MM 0.23 
C30-| M 0.14 


0.00 1.00 2.00 3.00 4.00 5.00 


Figure 2. Content of org. C, (%) in the soil profile 


It was found that in none of these layers the amount of org. C did not exceed 0.5%, 
which indicated that it could not be define as a spodic diagnostic horizon. The soil was 
determined as Cambisols. Nevertheless, there were processes, which indicated its future 
evolution towards podzolization, such as formation of acidic humus, acidification of the soil, 
destruction of clay minerals and transport of products in soil depth. 

A high amount of organic matter was accumulated in the surface soil horizon. The 
organic matter exceeded 15% in the first part of it (AA horizon). It should be noted that in the 
soil depth of the profile its amount decreased almost twice in the lower part of the A horizon, 
followed by a rapid decrease to EF and Bw horizons. This change was due to the removal of 
organic matter from EF horizon in depth. As a reason, high acidity can be noted. Soil response 
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is an important diagnostic indicator. In the Ah horizon it was estimated to be very highly 
acidic. The soil reaction had the lowest values in it (pHy20 = 4.1). The origin of this acidity 
was a result of the introduction of free organic acids from the decomposition of the roots that 
form the AA horizon and partly from the litter. According to some authors (Ganev, 1990) 
values of pHy20 below 4.8 are an indicator of the presence of free fulvic acids in the soil. 
Funakawa et al. (2006) point out that the pH around 4.0 is the critical limit that enhances the 
mobilization of org. C in the litter and in the surface soil layer. This causes its mobility into 
the soil depth of the profile in the form of a soluble organic substance. 

In the second part of the A horizon the value of pHy20 continued to be very low. This 
means that there is a removal of soluble organic matter from it. In depth of the soil profile, 
the pHy20 values increased slightly, in E horizon its values remained within the range of the 
highly acidic spectrum. The profile distribution of pHy20 showed the predominance of acidic 
products from the decomposition of the organic matter over the basic ones - of the weathering 
processes in the C horizon. 

The exchangeable acidity values, measured as pHcacn, showed that the soil buffer 
capacity was achieved due to Fe oxides dissolved as well as the weathering of secondary clay 
minerals (Ulrich scale, 1983). This referred to the depth from Ah to E horizon including. In 
the horizon, a major buffer process was the exchange of protons with basic cations. The 
titrimetric determined exchangeable acidity was very high. Its profile distribution is 
represented in Figure 3 and shows the high degree of acidification in the surface soil horizon, 
which is rich in humus. 


Exchangeable acidity, (cmol(+).kg-1) 


A res 
Aras 
E1420 
pw 20.20 
C30-) | 9 
Oo 2 4 6 8 10 12 


Figure 3. Exchangeable acidity in soil profile 


Under the influence of high acidity, the soil was highly leached. Basic cation values 
were low, with their lowest levels in the A horizon. The sum of the basic cations in the Ah 
and A horizon was: 1.99 cmol(+).kg™ and 0.42 cmol(+).kg’. Cation exchange capacity (CEC) 
values were highly differentiated between A, E and Bw horizons. In E horizon the CEC 
decreased 2.7 times compared to the A horizon. The results of the soil texture supported this 
assumption (Figure 4). 
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Table 2. Soil solution reaction, exchangeable cations, cation exchange capacity (T) and base 
saturation (V) 






































Soil pH pH IExch. |Exch. [Exch. | Exe [Exch. (Exch. | Exe | Sum of | CEC | BS 
Horizo | Depth acidity |Ca Mg h.K Na Fe h.M | basic T Vv 
ns n cations 

em  |H,0 (CaCl, [emol(+).kg! % 
An 0-7 4.1 B.3 16.80 1.29 0.45 0.23 0.02 0.27 0.01 | 1.99 19 11 
A 7-14 4.3 3.6 10.33 0.17 (0.13 0.11 0.01 (0.13 0.00 | 0.42 11 4 
E 14-20 4.7 4.1 4.05 0.17 0.07 0.04 0.02 0.03 0.00 | 0.30 4 7 
By 20-30 {5.2 4.3 4.76 7.00 {1.35 0.03 10.06 0.02 0.01 | 8.45 13 64 
Cc 30-| {5.4 4.3 1.99 6.22 1.59 0.03 10.07 0.02 0.01 | 7.90 10 80 


















































The profile distribution of the amount of clay showed its rapid decrease from A to E 
horizon - 6.7 times. The differences were slightly expressed in the soil depth. Bw horizon 
contained 0.6% less clay than EF horizon, at the expense of the amount of silt that increased 
1.4 times. The amount of sand in the £ horizon was commensurate with that in the C horizon. 


Sand, silt and clay distribution, (%) in soil 
profile 


SSS == a 55 370 
Ah 0-7 ARE 
_———SSSESEE eS fa 


A 7-14 — 


E1420 7.29 
ae za 
ai oes 
0 10 20 30 40 50 60 70 80 


sand silt m clay 


Figure 4. Sand, silt and clay distribution, (%) in the soil profile 


Thus, an F horizon was formed, which was clay-poor with an increased amount of 
sand. It had a coarse texture (sandy loam) than that in the Ah horizon (sandy clay loam). It 
could be assumed that the change of the texture began in the A horizon, where due to the loss 
of clay it was considered as sandy loam. The changes in clay content in the soil profile were 
accompanied by changes in cation exchange capacity. In E horizon CEC decreased rapidly - 
2.7 times compared to A horizon and 3.2 times compared to Bw horizon. This was due to the 
loss of colloidal particles - organic (Figure 2) and mineral (Figure 4). A podzolization process 
occurred in the soil. 

Some features of podzolization process were also shown by the Fe distribution. It was 
found to accumulate in Bw horizon, where its values reached 39 144 mg/kg (Figure 5). 
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Content of Fe, (mg/kg) 


Abo7 ea 

A714 

E1420 aa 

Bw 20-30 Ty 
C30-| ET ye 
Oo 5000 10000 15000 20000 25000 30000 35000 40000 45000 
Figure 5. Content of Fe, (mg/kg) in the soil profile 
Conclusion 


The analysis and assessment of the diagnostic features of the investigated soil did not 
prove the presence of a spodic diagnostic horizon, which did not allow the definition of 
Podzols. Processes such as: accumulation of acidic humus insitu, removal of org. C, 
acidification and loss of silicate clay, degradation clearly expressed by the accumulation of Fe 
in Bw horizon, were established. The evolution of Cambisols is related to a complex of 
specific processes, whose direction is pointed towards podzolization. 
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Abstract: 


The assessment of the distribution of soil physico-chemical properties facilitates the use 
of proper soil management practices. Hence, the study was conducted in Izombe rainforest of 
Nigeria to ascertain the profile distribution of physical and chemical soil properties. Samples 
were collected from two profile pits base on horizon differentiation for laboratory analyses using 
standard procedures. Data generated were analyzed statistically using coefficient of variation and 
correlation. The result indicated that sand recorded 760 g kg” in pedon 1, 730.20 g kg" in pedon 
2 and was predominant in the soil with higher amount at the surface horizon while clay and silt 
increased down the pedons. Bulk density ranges from 1.14 — 1.79 g cm®. The AB horizon was 
the least acidic of the pedons. Most soil chemical properties decreased down the profile. Sand 
and pH had low variation (4.80 % <12.60 %) while silt, organic carbon, total nitrogen had high 
variation (>41.52 % <79.58 %) among the horizons. Fe, Mn and Zn had moderate to high 
variation among the horizons. Clay correlated negatively and highly significantly (r= -0.829, r= - 
0.883, p= 0.01) available phosphorus and Fe. Clay also had a positive significant correlation (r= 
0.665, r= 0.641, p= 0.05) with bulk density and base saturation. pH had a significant negative 
relationship (r= -0.708, p= 0.05) with base saturation while cation exchange capacity had a 
highly significant positive correlation (r= 0.877, p= 0.01) with Pb. Fe had a negative and highly 
significant relationship (r= -0.838, p= 0.01) with Zn. The result indicated that most mineral 
element are domiciled in the A horizon and such demand proper management practices for soil 
nutrients sustainability. 


Keywords: distribution, horizon, rainforest, profile pit, soil properties, variation. 


Introduction 

A vertical cross section through the pedosphere shows varying degree of changes taking 
place at both surface horizon and subsurface horizons. These changes influence the utilization of 
these soils for agronomic, engineering, environmental and aesthetic uses. Variation in soil 
horizons exists at many scales with different dominant controlling factors. Parent material, 
climate, and geological history are major factors affecting the distribution of soil properties at 
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continental scale (Akamigbo and Asadu, 1983). Understanding the vertical distribution of soil 
properties along profile depths could help to predict the proneness of such soils to degradation. 
Knowledge of their variability is essential in applying location-specific land management 
interventions (Onweremadu ef al., 2015). Soil properties differ among soil profile horizons due 
to transport and storage of water and nutrient across and within the soil profile. Researchers (Oku 
et al., 2010) reported that soil pH increased with depth and it is least variable, irrespective of 
depth. They also found variability of organic matter, total nitrogen, available phosphorus, and 
CEC to increase with depth and ranged between moderate and high. 

Soil is likely to show great variability in their physical, chemical, biological and 
mineralogical properties because the soil is a heterogeneous unit. Knowledge of variability of 
soil properties is very indispensable as this determine the productivity and usage of the area. A 
study of the variability trends of soils is essential in order to highlight the soil potentials and 
enhance their management and productivity (Mahdi and Mark, 2006; Arnold, 1996). Ogunkunle 
(1993), working on the Alfisols of southwestern Nigeria, observed that soil pH was the least 
variable (low variability) property, irrespective of depth. Accordingly, horizons may differ in 
organic matter content, structure, texture, pH, base saturation, cation exchange capacity (CEC) 
and bulk density as well as many other soil physical and chemical properties. The rate of soil 
variation brings about rate of degradation from one location to another due to the composition of 
soils of the area. There is need to maintain the soil ecosystem through efficient land/soil 
management practices that will promote food security, soil nutrient balance and sustainable 
environment. Hence, the work is aimed to assess horizon variations in selected physical and 
chemical properties of Izombe soils. 
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Fig. 1. Location map of the study area (Google Imagery, 2018) 
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Materials and Methods 

Study area 

The study area as shown in fig | is IZombe in Oguta Local Government Area in Imo state 
south-east, Nigeria. It lies between latitude 5° 12’ and 5° 56’ N and longitude 6° 38’ and 7° 25’ E. 
The annual rainfall ranges from 2500 mm to 3000 mm with mean temperature of 28 °C (NIMET, 
2015).The original vegetation of the study area is the modified tropical rainforest which has been 
altered by human activities. The greater part of the area now is conserved with trees, shrubs, 
legumes. The geology of the area consists of alluvium which is marine deltaic deposits (Lekwa, 
1985). 

Field study and sample collection 

A reconnaissance study was carried out on the study areas. Subsequently, free survey 
technique was used to sink two (2) soil profile pits in the study area. Samples were collected 
following the guideline of FAO, (2006). Also, core samples were collected for bulk density 
determination. The samples were prepared for laboratory analyses. 

Laboratory analyses 

Particle size distribution was determined by hydrometer method (Gee and Or, 2002). 
Bulk density was determined by core method (Grossman and Reinsch, 2002). Saturated hydraulic 
conductivity was determined using a constant head method (Klute, 1986). Soil pH was determined 
using 1:2.5 soil — liquid ratio using a pH meter (Thomas, 1996). Organic carbon was determined 
by wet digestion method (Nelson and Sommers, 1982). Total nitrogen was determined by micro- 
Kjeldah digestion technique (Bremner, 1996). Available phosphorous was determined using 
Bray II method (Olsen and Sommers, 1982). Exchangeable acidity was determined by the 
method described by MacLean (1982). Exchangeable bases were determined by neutral 
ammonium acetate procedure buffered at pH 7.0 (Thomas, 1982). Available micro nutrients (Fe, 
Pb, Mn and Zn) were extracted by DTPA as described by Sahlemedhin and Taye, (2000) and all 
these micro nutrients were measured by atomic absorption spectrophotometer. 

Statistical analyses 

Coefficient of variation as used by Wilding et al. (1994) was used to estimate the degree 
of variability existing among soil properties in the study sites where CV< 15 % = low variation, 
CV =15 % <35 % = moderate variation, CV > 35 % = high variation. Correlation was also used 
to compare relationship between selected soil parameters. 


Results and Discussion 

The result (Table 1) indicated that sand fraction recorded (712 — 792) g kg"! and (682 — 
793) g kg in pedon 1 and pedon 2, respectively. Sand fraction also recorded low variation (4.80 
%, 9.20 %) in both pedons. The sand decreased down the profile with no specific trend in pedon 
1 while in pedon 2 it decreased in a specific trend. It was observed that the highest sand fraction 
was indicated in the Ap and AB horizons of pedon 1 and Ap horizon of pedon 2. The 
predominance of sand over other soil fractions could be attributed to parent material, climate and 
land-uses. This is in concurrence with the findings of (Obasi et al., 2015) on soils underlain by 
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coastal plain sand in some agro-ecology in southeastern Nigeria. The low variation could be 
associated with homogeneity of the soil. Soil with high sand content are associated with nutrient 
leaching, high infiltration and infertility. 


Table 1. Soil physical properties of the studied site 








Horizon Depth Sand Silt Clay TC Ksat BD Po 
cm g kg! cm/m gem” % 
Pedon 1 
Ap 0-30 792.0 40.0 168.0 SL 0.0334 1.14 54.81 
AB 30-60 792.0 20.0 188.0 SL 0.0198 1.52 59.60 
Bt 60-90 772.0 20.0 208.0 SCL 0.0114 ~=—1.64 48.25 
Bty 90-120 712.0 60.0 228.0 SCL 0.0194 1.58 55.18 
Bet 120-180 732.0 40.0 228.0 SCL 0.0430  ~=—-1.79 58.46 
Mean 760.0 36.0 240.0 0.0300 =1.53 55.26 
CV 4.80 46.48 12.78 49.69 15.78 8.01 
Pedon 2 
Ap 0-20 793.0 40.0 167.0 SL 0.0491 1.59 42.59 
AB 20-50 782.0 30.0 188.0 SL 0.0480 1.59 42.40 
Bt 50-75 752.0 20.0 228.0 SCL 0.0582 1.69 21.86 
Btp 75-105 692.0 60.0 248.0 SCL 0.0590 ~—:1.70 21.86 
Bet 105-180 632.0 10.0 268.0 SCL 0.0690 —1.70 21.56 
Mean 730.2 32.0 219.8 0.0600 =1.65 30.05 
CV 9.20 60.11 19.01 15.09 3.54 37.79 





TC= textural class, SL= sandy loam, SCL= sandy clay loam, Ksat= saturated hydraulic conductivity, BD= 
bulk density, Po= porosity 


Silt fraction recorded high variation in both pedons (46.48 %, 60.11 %) but also had no 
specific trend of decrease with soil depth. The high variation could be associated to intense 
rainfall that facilitated sorting of soil particles. The silt content is in conformity with the findings 
of Madueke et al. (2012) that soils of southeastern Nigeria are low in silt as a result of the high 
degree and extent of weathering and leaching they have undergone. 
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Clay fraction (Table 1) recorded (168 — 223) g kg” and (167 - 268) g kg" in pedon 1 and 
2, respectively. However, clay increased down the profile while it recorded low variation in 
pedon 1 and moderate variation in pedon 2. The clay content of the soils increases with the 
increasing rainfall and annual temperature. The variation could be attributed to rate of clay 
migration, drainage pattern and soil structure. However, Obi and Akinbola (2009) reported that 
increase in clay with soil depth could be attributed to clay translocation and erosion in surface 
horizon. This is similar to the findings of (Osujieke et al., 2017) in soils of southeastern Nigeria. 
It was observed that the Bt horizon had high clay content which is an indication of its argillic 
properties. Higher clay content in the subsurface than the surface horizons in the pedons with 
subsequent formation of argillic subsurface horizon may be attributed to sorting of soil materials 
by eluviations and illuviation processes of soil formation. 

Kya recorded mean of 0.03 em m” and 0.06 cm m’ for pedon 1 and 2 but was very slow 
according to the rating (<0.2 cm m’) of Landon, (1991). Ksat is influenced by soil texture and 
structure which in turn is influenced by total porosity and pore size distribution. Also, Kgat 
recorded high variation (49.69 %) in pedon 1 and moderate variation (15.09 %) in pedon 2. The 
rate of variation could be associated with soil texture, bulk density, available pore spaces and 
intense rainfall. The result showed that soils of Bct horizon conducts more water over other 
horizons in both pedons. The implication of slow hydraulic conductivity in soils is that the soil 
will have high runoff and more erosion. 

Bulk density (Table 1) recorded (1.14 - 1.79) g em® and (1.59 - 1.70) g cm® in pedon 

land 2, respectively. Also, bulk density increased down the profile but with an irregular trend in 
pedon 1. The bulk density values in this soil unit were lower than 1.85 g cm”, a value considered 
as very high by SSS, (2006) to impede root penetration. However, the value of bulk density is 
low when compared with findings of (Onweremadu et al., 2011) in soils of southeastern Nigeria. 
However, the lower values at the surface soil were probably due to presence of organic matter 
content. The increase in bulk density with profile depth could be associated to decrease in 
organic matter content, less aggregation and compaction caused by the weight of the overlying 
soil layers, presence of parent material and low faunal and root activities. Low bulk density 
influences high porosity which is an indication that soils are not compacted. 
Soil pH (Table 2) was slightly-strongly acidic in both pedons according to the ratings of Chude 
et al. (2011). The pH had an irregular pattern of distribution with soil depth. The AB and Bet 
horizon is the least and most acidic, respectively. Soil pH recorded low variation (12.60 %, 9.52 
%) in both pedons which may have resulted from the parent material and intense rainfall which 
wash off most basic cations. Several researchers (Ndukwu et al., 2012) have reported on low 
variation of soils of eastern Nigeria and it is the consequence of acidic nature of the parent rocks, 
coupled with the influence of the leached profile under high annual rainfall condition. 

Organic carbon was low according to the ratings of (Akinrinde and Obigbesan, 2000) and 
recorded range of 0.14 % - 0.64 % and 0.18 % - 0.63 % in pedon | and 2. However, the level of 
organic carbon of the studied pedons was low compared with the findings (Amanze et al., 2016) 
on soils of southeast Nigeria. Also, Ogeh and Ukodo, (2012) reported low OC in soils of 
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rainforest zone of Nigeria. Organic carbon (Table 2) recorded high variation (56.42 %, 79.20 %) 
which could be attributed to higher organic matter deposit on soil surface, rate of decomposition 
and leaching due to intense rainfall. This is conformity with the findings of (Ukaegbu et al., 
2015) on soils of southeast Nigeria. However, OC decreased down the profile but with an 
irregular pattern in pedon 2. The Ap and By horizon recorded the highest level of OC in pedon 1 
and 2, respectively. The OC content on Ap horizon could be attributed to the organic material 
deposited on the surface horizon while that of the By: could be as a result of other carbon 
attributors in the soil. 

Total nitrogen had range of 0.056 % - 0.07 % and 0.028 % - 0.098 % in pedon 1 and 2, 
respectively. The total nitrogen level was low according to the ratings of (Enwezor, 1990). The 
low content of total N in the soils could be attributed to low organic matter of these soils, since 
inorganic N is accounting for only a small portion of total N in soils (Almu and Audu 2001). The 
low amount of total N reflects the amount of organic carbon in the soils. Nitrogen can be 
introduced into the soils by natural processes such as lightening, decay of animal and plant 
tissues while the main cause of N deficiency in tropical soils is intense leaching and erosion due 
to the intense tropical rainfall. 

Available phosphorus recorded mean of 21.74 mg kg" and 18.25 mg kg" in the studied 
pedons and were high according to the ratings (>15 mg kg’') of (Enwezor, 1990). Available P 
was high compare with the findings of (Osujieke et al., 2016) in soils of southeast Nigeria. The 
available phosphorus recorded moderate variation (20.70 %, 23.09 %) in both pedons which 
could be attributed to organic material deposit, rate of fixation and leaching. The result is 
contrary to the findings of (Busari et al., 2005; Uzoho and Oti, 2005) that have reported high P 
deficiency for tropical soils. 

Exchangeable bases (Ca, Mg, Na, K) recorded mean of (2.25, 1.29, 0.095, 0.104) cmol 
kg" in pedon 1 and (2.64, 1.40, 0.08, 0.09) cmol kg" in pedon 2. Ca was low, Mg was moderate, 
Na and K were very low according to the rating of Landon (1991). The levels of exchangeable 
cations are dependent on the nature of parent materials and extent of weathering. The 
exchangeable bases were within the range reported by Amanze et al. (2016) in soils of southeast 
Nigeria. Exchangeable cations had no specific trend of distribution down the profile pits. The Ca 
and Mg had low variation in both pedons while Na and K had low variation in pedon 1 and high 
variation in pedon 2. The variation among the exchangeable cations is associated with leaching 
and runoff wash of the basic cations, pH level of the soil and fixation rate of the cations. 

Cation exchange capacity had range of (2.20 — 4.40) cmol kg” in pedon 1 and (3.00 — 4.40) cmol 
kg" in pedon 2. CEC was low according to the ratings of (Esu, 1991) on tropical soil. The CEC 
result was low compared with the works of (Amanze et al., 2016) on soils of southeast Nigeria. 
The low CEC could be due to the domination of low activity clay such as kaolinite in these soils 
which resulted from highly weathered parent rock. This is in concurrence with the finding of 
Akinrinde and Obigbesan, (2000) on soils of southern Nigeria. Also, Ukpong, (2000) stated that 
low level of cation exchange capacity in soils could be associated with tidal imports, runoff and 
seepage. However, CEC decreased down the profile in no specific trend. The Ap horizon 


Bulgarian Journal of Soil Science® 2018 Volume 3. Issue 2. www.bsss.bg 95 


Osujieke et al., 2018 





recorded the highest level of CEC which suggested that organic matter could have been its major 
attribute. The CEC had moderate variation (29.06 %, 16.40 %) in pedon | and 2 which could be 
associated with organic matter level and runoff. Hence, (Enwezor, 1990; Esu, 1991) reported that 
organic matter is the major attribute of CEC in soils of southeastern Nigeria. Lombin et al. 
(1991) also reported that organic matter content was a major contributor to the CEC of the soil. 
Therefore, there is need to incorporate organic material into the soil to enhance the exchange 
site. 

Available Fe was below the critical limit recommended by (Kparmwang et al., 2000, Deb 
and Sakal, 2002). The Fe content was low compared with the findings of (Ahukaemere et al., 
2016; Ndukwu et al., 2016) on soils of southeastern Nigeria but was high when compared with 
the findings of (Mulima et al., 2015) in soils of Northeast Nigeria. The coefficient of variation 
showed that available Fe recorded high variation (37.82 %) in pedon 1 and moderate variation 
(32.77 %) in pedon 2. The variation could be associated with rate of mineralization and intense 
rainfall activities which is associated with study area. Fe content decreased down the profile but 
had no specific trend in pedon 1. 

Lead (Pb) recorded (0.089 — 0.630) mg kg" in pedon | and (0.216 — 0.658) mg kg" in 
pedon 2. The Pb content was above the critical limit recommended by (Deb and Sakal, 2002; 
Esu, 1991). Pb recorded high variation (72.67 %, 46.57 %) in both pedons which be attributed to 
the fixation rate and the amount released by organic matter and mineral matter. The Pb decreased 
down the profile in an irregular trend. The Pb level was lower than the result reported by 
(Osakwe and Okolie, 2015; Das et al., 2015) in soils of southern Nigeria. The lead level of the 
soil is not at toxicity level and thereby will not be harmful to the soil 
Mn recorded mean of 2.56 mg kg”! and 2.27 mg kg” for soils under pedon 1 and 2, respectively. 
Mn was below the critical limit recommended by (Sims and Johnson, 1991). Hence, the available 
Mn is suitable for crop production and environmentally friendly. It had an irregular pattern of 
decrease with soil depth and recorded the highest level in Bt] horizon in Pedon | and Ap horizon 
in pedon 2. These values are within the range values reported by (Udoh et al., 2008) in soils of 
southern Nigeria. However, Mn had moderate variation (27.53 %) in pedon | and low variation 
(10.86 %) in pedon 2. The variation could be attributed to soil textural class and rate of porosity. 
The AB horizon had the least Mn in both pedons while the By; horizon recorded the highest level 
of Mn in pedon 1 while the Ap horizon recorded the highest Mn level in pedon 2. 

Zinc recorded range of (0.245 — 0.640) mg kg” in pedon 1 and (0.247 — 0.520) mg kg" in 
pedon 2 with no specific trend of increase down the profile. Zn level is below the critical limit 
according to the ratings of (Esu, 1991; Kparmwang et al., 2000; Deb and Sakal, 2002) which is 
an indication of non toxicity in the soil and as such will pose no threat to the environment. The 
value of Zn are below the range reported by (Udoh et al., 2008) in soil of southeastern Nigeria 
but above the range reported by Oyinlola and Chude, (2010) in soils of Northeastern Nigeria. 
The available Zn was more predominant at the B horizons which could be attributed to 
illuviation of materials. 
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The sand fraction had positive and significant correlation (t= 0.753, p= 0.05) with the 
available phosphorus while it had negative and highly significant relationship (r= - 0.934, p= 
0.01) with the clay fraction. The clay fraction had positive and significant relationship with the 
bulk density (r= 0.665, p= 0.05) and the base saturation (r=0.641, p= 0.05), while it had negative 
and highly significant relationship with available phosphorus (r= - 0,829, p= 0.01) and Fe (r= - 
0.883, p= 0.01). Ksat had negative and highly significant relationship (r= -0.804, p= -0.01) with 
total porosity. This relationship may be associated to the textural composition of the soil. 

However, Brady and Weil (2002) stated that the rate of hydraulic conductivity is 
dependent on soil pores. The pH had negative and significant relationship (r= -0.708, p= 0.05) 
with base saturation. This relationship could be associated to the leaching of most of the base 
forming cations thereby leaving the exchange complex dominated by acidic cations. However, 
the pH had negative relationship (r= -0.234, r= -0.230, r= -0.267) with Pb, Mn and Zn but related 
positively (r= 0.179) with Fe. Organic carbon correlated negatively and significantly (r= -0.713, p= 0.05) 
with total nitrogen. This may be due to the factor of age of the soil organic matter, the prevailing climatic 
condition, N-leaching and volatilization and microbial N-assimilation. The CEC correlated positively and 
highly significantly (r= 0.877, p= 0.01) with Pb. Bulk density correlated significantly positive (r= 0.665, 
p= 0.05) with clay but had negative correlation (r= -0.117) with OC. However, Fe had negative and 
highly significant relationship (r= -0.883, r= -0.747, r= -0.838, p= 0.01) with clay, base saturation and Zn. 
However, the positive correlation implies that increase in one soil property increases the other 
while, negative correlation implies that increase in one soil property decreases the other and vice 
versa 

Conclusion 

This study has shown soil physicochemical properties and the variation among the 
horizons of the studied pedons. The organic carbon, total nitrogen, and exchangeable basic cation 
(Ca, Na, K) were predominantly low while available phosphorus was high. Also, Fe, Mn and Zn 
were below the critical limit required in the soil. Soil properties most be considered for 
sustainable agriculture and proper land management practices that will promote the environment. 
Application of organic matter will be essential for most of the nitrogen and basic cation to be 
utilized by crops for quality productivity. 
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Abstract 

This field study shows how different crops and soil treatments with different nitrogen 
fertilization levels affect sulphur balance. Sulphur deposition, leaching through soil and water 
seepage, availability and uptake by crops, effects on crop yield and losses were investigated. 
The research was conducted on a field trial from 1996 to 2014 in temperate continental 
climate on Stagnosols. The soil was treated 10 times with mineral nitrogen fertilizers in an 
ascending doses from 0 to 300 kg N ha’ and additional two 250 kg N ha’ treatments, one 
with phosphogypsum and one with zeolite tuff and CaCO3. Drainpipes and lysimeters were 
installed to collect water samples and measure the effects of N doses, amendments and 
precipitation on sulphur losses. Average annual sulphur losses were from 4.9 to 68.7 kg ha 1 
through drainage water, and from 1.5 to 24.9 kg ha’! through lysimeter water. Depending on a 
crop, year and yield, average crop losses ranged from minimal 2.8 kg ha” in winter wheat up 
to 17.6 kg ha’' in oilseed rape. Average total S content in soil varied from 882 to 1764 kg ha’. 
Overall agroecosystem S balance between input and output, calculated for all crops and all 
treatments was positive only in the treatment with phosphogypsum for winter wheat, corn and 
oilseed rape. Water shortage in 2011 caused positive S balance for soybean in other 
treatments too, not only in the treatment where phosphogypsum was applied. 


Key words: phosphogypsum application, mineral nitrogen fertilization, crops, 
drainage and lysimeter water, soil, atmospheric deposition 


Introduction 

The presence of sulphur in the soil is in the form of organic compounds, sulphides 
(S~), elemental sulphur (S°) and sulphate (SO,7-). According to Isaac, Kerber (1971), a 
normal concentration range of total sulphur in soils is from 500 mg kg" to 4000 mg kg". 
Total sulphur in mineral soils may range from < 20 mg kg” in sandy soils up to > 600 mg kg 
in heavy texture soils. Organic soils may contain as much as 5000 mg kg” (Zgorelec et al., 
2012). Most soils, however, contain between 100 and 500 mg kg’! of sulphur. Elemental 
sulphur must be oxidized by soil bacteria to SO,7~ before becoming plant available. Warm 
temperatures and good moisture and aeration are required for S-oxidizing bacteria to work. 
Sulphur oxidation is minimal at soil temperatures that are less than 10 °C. Even at 23 °C, 
sulphur oxidation rate is about 15 % of the oxidation taking place at 29 °C, therefore the peak 
rates of sulphur oxidation do not occur until late spring (Camberato et al., 2012). After 
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sulphur is oxidized into sulphate, plants absorb it and it is converted into amino acids cysteine 
and methionine. These amino acids are the building blocks of sulphur based proteins and 
other key role compounds in multiple plant functions like photosynthesis, disease resistance 
and plant growth. Because of its important roles, optimizing the sulphur nutrition of crops is 
key to achieve top crop yields and quality. Since the early 1960s, the research on the response 
of crops to sulphur expanded over years and covered a large number of crops, soils and 
fertilizers. By 1990, it was clear that sulphur deficiencies were becoming increasingly 
important as evidenced not only by soil analyses but, more importantly, by the results of field 
research which documented effects of sulphur application to about 30 crops under various 
weather conditions and on various types of soil (Tandon, 2011). Sulphate leaching may be 
subject to considerable variations caused by differences of pH in the soil. Soils with lower pH 
value have better adsorption of sulphur. Eriksen et al. (2002) investigated sulphate leaching 
and S balance in organic cereal crop rotation. It was found that sulphate leaching was 
quantitatively the most important item of the S balance. On some treatments of nitrogen 
fertilization, FG was used. In the process of decomposition of raw phosphate with sulfuric 
acid which was used for phosphoric acid production, FG occurs as by-product. Vyshpolsky et 
al. (2010) conducted a two-year experiment on the plots of cotton on heavy clay with a 1% 
content of organic matter, and by applying FG the efficiency of irrigation was increased and 
water retention around the root system improved. 

Given the fact that sulphur has been widely recognized as the fourth major plant 
nutrient after nitrogen, phosphorus and potassium and cannot do the job alone than it must 
work as a member of the plant nutrition team, assessment of the sulphur balance in this 
experiment is significant and gives an idea of the sulphur cycle in the soil-fertilizer-plant 
system. Sulphur addition with phosphogypsum, amelioration with zeolite tuff and CaCO; and 
fertilization with different nitrogen levels allows the role of sulphur and nitrogen to be 
monitored for growth, yield and crop quality, to be assessed the impact of pH on sulphur 
mobility and extraction, and of the plants, as well as sulphur adsorption from the soil. On the 
other hand, an anthropogenic impact through human influence on agriculture affects the 
ecosystem in general. Excessive accumulation of sulphur in the soil represents an 
environmental problem through its deposition on the soil and it’s leaching in the lower layers 
of the soil which can lead to soil contamination. Therefore, it is important for us to explore the 
ecosystem, to have the agroecosystem and to implement management of plant nutrition, but in 
the same time to protect the ecosystem. Accordingly, the main goal was to determine sulphate 
losses and gains in a field experiment with a certain soil type and climate including different 
nitrogen fertilization doses and to valorise and determine an overall sulphur balance in 
agroecosystem (soil, planted crops, water and air deposition). 


Materials and Methods 

This research was conducted from 1996 to 2014 in central part of Croatia in Popovaca 
(Figure 1) where intensive agriculture occurred. Selected location is on the border with Park 
of nature Lonjsko polje which has natural water retention and needs to be protected. 
According to Skoric (1986), the type of soil was plain, deep, district drained pseudogley 
(Stagnosols); deep, district pseudogley according to Husnjak (2014), and district Stagnosols 
according to FAO, WRB (2006). Excessive humidity was present in the upper part of profile, 
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because of rainfall, and ground water occurred below 175 cm of depth. The soil was slightly 
porous throughout its depth, with mediocre capacity for water. Air capacity was very small. 
The soil was classified as acid (~ 4.5) (Zgorelec et al., 2013). The compaction of soil in the 
field trial was medium in arable and eluvial, and higher in deeper horizons. Because of these 
characteristics, drainage pipes were conducted on the surface. The characteristic of the area is 
continental climate with mean annual temperature of 10.7 °C and average rainfall of 865 mm 
(Zgorelec et al., 2013). 
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Figure 1. Study location 


The trial consisted of 10 treatments with different mineral nitrogen doses. Table 1 
presents treatments with different fertilization doses and application time during the 


experiment. 


Table 1. Fertilization (kg ha’) treatments and application time (per year) 





























Treatment Fertilization (kg ha") Application period/year 
I. Control - no fertilization - 
Il. ON+120P+ 180K every year 
Ill. 100N+120P+ 180K every year 
IV. 150N+120P+ 180K every year 
V. 200 N+ 120P+ 180K every year 
VI. 250N+120P+ 180K every year 
-1 
VIL 250N+120P+180K+12tha 1996, 2002 and 2009 
Phosphogypsum (FG) 
VIIL. 250N+120P+180K+3 tha" Zeolite 1996, 1999, 2002, 2009, 
tuff + CaCO; (+Dolomite) (+2013) 
IX. 300 N+ 120P+ 180K every year 
Xx. Black fallow - 





All agro technical measures including tillage, herbicide applications etc. were identical 
in all trial treatments. 

Soil samples were taken three times per year (2010), at the beginning of the growing 
season, during the growing season and immediately after the harvest of certain crops (winter 
wheat). Samples were taken from each treatment at surface horizon (0 — 30 cm). Total S in 
soil (kg ha”') was calculated from the measured concentration (mg kg”) of S in soil and soil 
bulk density (1.4 g cm™) determined in arable horizon (0-30 cm). The measurement for total S 
content in plant material started in 2006. The losses of total sulphur in kg ha’! were calculated 
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from measured concentration of total sulphur (mg kg") in plants and from measured crop 
yield in kg ha’'. Lysimeters and drainpipes were installed as described in Zgorelec et al. 
(2013); Vukovic et al. (2008) and Sestak et al. (2014). In drainage and lysimeter water, S- 
SO, concentration in water (mg L”) and drainage discharge (L day’) are taken to calculate 
sulphur losses in kg ha’ (Figure 2 and Figure 3). Table 2 presents parameters determined in 
soil, plant material and water including determination methods. 


Table 2. Parameters determined in soil, plant material and water including determination 
methods 


























. Soil quality - Sampling - Guidance on HRN ISO 
Sampling, transport and . ; 
cosas sampling techniques 10381-1 to 8 
S (2002-2006) 
. er Soil samples preparation for physical HRN ISO 
d /grinding/ 
: hie eo as ‘ and chemical analysis 11464:2006 
sieving/homogenizing 
SOIL 
Soil quality - Determination of pH 
HRN I 
pH in 0.01 M CaCl, 1 M KCl and H,0 at S 
; 10390:2004 
1:2.5 (w/v) ratio 
TS Soil quality - Determination of total HRN ISO 
(% DM; g kg”) sulphur with dry combustion method 15178:2005 
sampling Cereals and cereal products - sampling ISO 24333:2009 
; se ; : ified HRN 
drying/grinding/ Sample preparation for analysis of haces a 4-2006 
PLANT sieving/homogenizing inorganic chemical composition 
TS Determination of total sulphur with dry = modified HRN 
(% DM; g kg”) combustion method ISO 15178:2005 
Water quality - Determination of 
fluorides, chlorides, nitrites 
F, Cl, NOy, Br, NO;, SO,” ; oes 
WATER a oO ; : orthophosphates, bromides, nitrates and HRN ISO 


sulphates with ion liquid 10304-1:1998 
chromatography - part 1. Method for 
water with low contamination 


(mg kg"') 





The crops grown in crop rotation are presented in Table 3. Considering that 
measurement of total S content in plant material started in 2006, the collected, analysed and 
taken into balance account data are marked with the bold text in Table 3. 


Table 3. Cultivated crops during the experiment 








Crop Vegetation year 
Corn (Zea Mays L.) 1996, 1999, 2004, 2007, 2013 
Oilseed rape (Brassica napus L.) 1998, 2001, 2009 
Soybean (Glycine L. max) 2002, 2005, 2011 
Winter wheat (Triticum aestivum L.) 1997, 2000, 2003, 2006, 2008, 2010, 2012 
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In this research, quality assessment and control were included. Analytical laboratory of 
Department of General Agronomy integrates Internal and External quality control. Internal 
quality control is carried out on a daily basis and involves checking the accuracy and 
repeatability of measurements by analysing the reference material (RM) and doing a 
minimum of three consecutive measurements. RM samples from different inter laboratory 
trials were used in the laboratory: ISE (International Soil analytical Exchange Programme 
organized by WEPAL [Wageningen Evaluating Programs for Analytical Laboratories, 
organized by Soil Chemistry and Chemical Soil Quality group, Department of Soil Quality of 
Wageningen University (WU), Netherlands]) for soil, IPE (International Plant analytical 
Exchange Programme organized by WEPAL) for plant and IFA (Water Proficiency Testing 
scheme organized by BOKU, Department for Agrobiotechnology, IFA-Tulln, Canter for 
Analytical Chemistry) for water. External quality control is performed once a year where it 
participates in various plant, soil and water field testing programs. 


Results and Discussion 

Drainage water 

Figure 2 shows minimum, maximum and average annual SO,-S losses through 
drainage water which are calculated for each treatment for the period between 1996 and 2014. 
If treatments VI, VU and VIII, to which the same doses of N, P and K were applied, are 
compared and if the fact that 335 kg ha’ of S was annually applied through FG in treatment 
VII, is taken into consideration, the average loss was the highest in treatment with FG (VID, 
which is 9.4 times higher than the average loss in treatment VI leads to the conclusion that 
losses of S due to FG application were 61.4 kg ha’'. Furthermore, we revealed that 18.3 % of 
applied S (in treatment VII) was drained, while the rest was adsorbed by the plant or the soil. 
Accordingly, Monaghan et al. (2000) investigated nutrient losses in drainage water where 
leaching losses of sulphate in the 0 N fertilizer treatment was 1.9 times higher (17 kg ha’') 
than in this research. In Riley et al. (2002) research, three years experiment recorded total S 
loss of 34.9 kg ha’. 
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Figure 2. Average annual SO4-S losses through drainage water (from 1996 to 2014); n=151 
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Lysimeter water 

Figure 3 shows minimum, maximum and average annual SO,-S losses through 
lysimeters water which were calculated for each treatment for the period from 1996 to 2014. 
If treatments VI, VII and VIII, to which the same doses of N, P and K were applied, are 
compared and if the fact that 335 kg ha” of S was applied annually through FG in treatment 
VII, is taken into consideration, the average loss was the highest in treatment with FG (VI), 
which is 9.7 times higher than the average loss in treatment VI meaning that the input of S to 
water due to FG was 22.4 kg ha’'. The conclusion is that 6.7 % of applied S gets lost through 
lysimeter water, while the rest is adsorbed in the soil and/or the plant. Eriksen, Askegaard 
(2000) found that sulphate leaching from an organic dairy crop rotation, on sandy soil was 20 
kg ha! — average of 4 years, equivalent to 60 % of total sulphur input. In this research, 
average S loss by lysimeter water calculated for all treatments except treatment VII was 2.0 
kg ha’. 
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Figure 3. Average annual SO4-S losses through lysimeter water (from 1996 to 2014); n=124 


Losses of sulphur through crop yields 

It is well established that sulphur application not only increases crop yield but also 
improves crop quality. It improves crop quality in many ways, some of which are: increasing 
the oil content of seeds, improving nutritional quality of forages by providing a balanced NS 
ratio, increasing crude protein content, and improving baking quality of wheat (Tandon, 
2011). A balanced NS ratio is important because of nitrogen use efficiency, plant vigour, 
water use efficiency, phosphate use, carbohydrate production and utilization, the rate of grain 
fill, maturity and many other plant factors. It primarily reflects the complementary 
relationship that nitrogen and sulphur have in producing plant proteins. On a higher nitrogen 
fertilization level, the plant uptake of sulphur was higher, because on a high nitrogen level in 
the soil, sulphur in the soil becomes more mobile, and the plant can easily take it (synergism). 
The importance of concurrent NS management in wheat was described in Salvagiotti et al. 
(2009) research. In this research, the highest level of nitrogen fertilizers was applied in 
treatment [IX with 300N + 120P + 180K kg ha’'. Losses through yield were calculated by 
sulphur content in crop yield. The lowest S losses were observed in winter wheat at treatment 
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with ON + 120P + 180K kg ha’ (II), and the highest in oilseed rape at treatment with 250N + 
120P + 180K kg ha’! + zeolite tuff + CaCO3 (VIII) (Figure 4). The sulphur requirements of 
different crops vary widely. It is reported by Beeson (1941) that legumes utilize large amounts 
of sulphur for growth, while cereals and similar crops generally have less requirements of 
sulphur. The author also revealed that sulphur content in wheat (grain/straw) ranged from 1.7 
g kg to as high as 4.5 g kg" in sweet clover. In this research, sulphur was also associated 
with a higher biological N fixation by legumes. The trend of increase in S content in grain due 
to rising mineral N content in the soil was noticed in Figure 4. For example, in soybean, the 
application of sulphur increased the number of nodules/plant and nodule weight (Gupta, 
Dubey 1998; Tiwari, Gupta 2006). Whenever the sulphur status of growing plants drops 
below the critical level which is different for different cultures, the visual symptoms of 
sulphur deficiency start appearing on the plant, generally in leaves. In many ways, sulphur 
deficiency symptoms resemble those of N in the early stages in that the leaves become pale- 
yellow or light green. The appearance of such symptoms itself indicates a serious condition 
because crop yields can decrease even without the appearance of such symptoms as a result of 
hidden hunger (Tandon, 2011). In Figure 4 losses of total sulphur through different crop yield 
are shown respectively (winter wheat, corn, soybean and oilseed rape). 

Table 4. Average sulphur content (g kg” ), annual average yield (t ha’) and losses of S 
through yield for all treatments (kg ha") 








Winter wheat Corn Soybean Oilseed rape 
n=144 n=72 n=36 n=36 
OYEIAEE salphe : 20 1.8 1.7 4.9 
content in yield, g kg 
Annual averse yield, t 3.2 6.7 2.7 2.6 
ha 
Average losses of S ea 12.1 4.6 12.7 


through yield, kg ha" 





In Table 4, average sulphur content in yield, annual average yield and losses of S 
through yield were calculated for all treatments and all years. Table 4 clearly shows higher S 
content in oilseed rape than in other crops, while losses of S in oilseed rape were similar to 
corn, because of higher yield of corn.Sulphur removal by corn, winter wheat, oilseed rape and 
soybean are presented on Figure 4. Average losses of S through winter wheat yield (Table 4) 
is in consistence with Beeson’s (1941) meaning that harvesting crops with lower percentage 
of sulphur content may result in annually removal of 6 kg ha’! of sulphur. Sulphur deficiency 
significantly effects the production and quality of winter wheat (McGrath, 2003; Gyéri, 
2005). Without adequate sulphur, crops cannot reach their full potential in terms of yield, 
quality or protein content; nor can they make efficient use of applied nitrogen (Sahota, 2006). 
Eriksen et al. (2002) reported that sulphur output by winter wheat uptake could be up to 4 kg 
ha’' depending on year and soil type. 


Bulgarian Journal of Soil Science® 2018 Volume 3. Issue 2. www.bsss.bg 110 


Galic et al., 2018 





According to this research, in Camberato, Casteel (2017) research about 11 kg S ha” is 
removed by corn. The research of sulphur application to corn and soybean crops in Iowa 
revealed that the combined S supply from soil profile, sulphate, atmospheric deposition, crop 
residue and organic matter mineralization fulfil corn requirements for sulphur (Sawyer et al., 
2002). 

For soybean, sulphur requirement is higher especially than in many other crops (Akter 
et al., 2013). Considering that fact, losses in this research are low compared with other 
investigated crops (corn, winter wheat and oilseed rape). Received losses can be attributed to 
yield, which was considerably lower in 2011 due to the intense drought. In Lamond (1997) 
research, sulphur removal by soybean grain was 13 kg ha”. 

In Eriksen, Mortensen (1999) research, removal of S by oilseed rape was 3.6 and 3.0 
kg ha’', respectively. Oilseed rape has a high demand for sulphur compared with cereals, with 
approximately 16 kg of sulphur required to produce each tone of seed (HGCA, 1998). 
Consequently, oilseed rape is particularly sensitive to any shortfall in sulphur supply. 
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Figure 4. Losses of sulphur through crops yield/grain (corn, n=8; winter wheat, n=16; 
oilseed rape, n=4; soybean, n=4) 

Sulphur content in soil 

The main source of S for plants originates from soil. According to Yang et al. (2007), 
Kotkova et al. (2008), and Forster et al. (2009), the total soil S content depends on fertilizer 
types applied (mineral or organic manure), and changes in soils have been found to be 
proportional to the amount of organic residues added. As it is written in the introduction, a 
normal concentration range of total S content in soils could range from 2 100 kg ha’ to 16 
800 kg ha™', while in this study the range was from 882 kg ha’ in control treatment (I), to 
1764 kg ha’ in 250 kg N ha’ on FG variant (VII) (Figure 5), in soil layer 0 — 30 cm. In total, 
36 000 kg FG ha’ was applied three times every six years, i.e. 12 000 kg FG ha’ per dose, on 
treatment VII. Accordingly, the S content that was applied in this study for a total period was 
6 698 kg S ha’, or 335 kg S ha” per year. If we compare treatments VI and VII, on which 
were applied the same doses of NPK but they differ in FG, we observed 462 kg ha’! more S in 
the soil on treatment VII with FG than on treatment VI without FG. Considering that FG 
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contains S, it can be added through the FG application where the exact concentration of S 
must be determined according to its leaching tendency. In a long-term field experiment 
conducted by Yang et al. (2007), total S ranged from 1392 kg ha’ in the control without 
fertilizer application to 1808 kg ha’ in the treatment with farmyard manure. In the plots 
receiving either NK or NPKS fertilizers for more than 80 years, no significant accumulation 
of total S, as compared with the control, could be detected. 
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Figure 5. Average sulphur content in soil 


Sulphur wet deposition through trial period 

Deposition of sulphur compounds is considered a major environmental issue 
(pollutant) — especially in Europe, where its effects were first observed (Walaszek, 2013). 
Due to technological development in industry and switching fuels from coal to oil and gas, 
emissions of oxidized sulphur dropped by 45 % over the last few decades. The removal of 
pollutants from the atmosphere can occur in one of two ways: by wet deposition or dry 
deposition (Fagerli, 2008; Giannitrapani, 2006). Erisman et al. (2005) estimated that about 80 
% of pollutants deposited in Europe are subject to wet deposition. Wet deposition is a process 
of washing out pollutants from the air by precipitation or cloud droplets (Walaszek, 2013) and 
may cause natural vegetation damage (Bytnerowicz et al., 2007), soil and water acidification 
(Sullivan et al., 2005; Zhao et al., 2009) and water eutrophication (Camargo, Alonso 2006). 
As it stands in Hu et al. (2005) current inputs of S from atmospheric deposition are less than 
10 kg ha’ in most Western European countries. Annual wet deposition of S— SO,” at 
measuring station Bilogora, which is the nearest to the field trial ranged from 2.6 kg S — 
SO,” ha’ in 2011 to 10.6 kg S—SO,? ha’ in 2001 (Statistical Yearbook of Republic of 
Croatia, 2001 - 2015) and shown on Figure 6. An average deposition through this study was 
6.0 S—SO,4* ha. Jian et al. (2014) studied atmospheric wet deposition of sulphur in 
developing and developed areas of South-eastern China and found S deposition in range of 
56.02 — 59.06 kg ha ' per year, respectively, surpassing their corresponding critical loads in 
China which is incomparably higher than in this research. Knights et al. (2000) reported that 
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large S inputs from atmospheric deposition and from sulfate fertilizer did not accumulate as 
organic and inorganic S in arable soils, but organic S obviously increased after FYM supply. 





12 


kg ha-' 














Figure 6. Annual wet deposition of S-SOq from 2001 to 2014 (Source: Statistical Yearbook of 
Republic of Croatia; 2001] - 2015) 


Sulphur balance through trial period 

Sulphur balance for all investigated crops is shown in Table 5. Volatilization from 
crops or soil was ignored as the amount of sulphur emitted from these sources was considered 
insignificant (Janzen, Ellert 1998). The input from atmospheric wet deposition was in average 
6.0 kg ha’! of S-SOq. The following year, an average sulphur balance was calculated for 
specific crops: winter wheat — 2006, 2008, 2010, 2012; corn — 2007, 2013; soybean — 2011; 
oilseed rape — 2009. The results of average losses of drainage water, lysimeter water and 
removal through yield as output were shown for every year that was taken in calculation. The 
average wet deposition and FG application were taken as input value of sulphur. The overall 
sulphur balance was negative, except in treatment where FG was applied. 
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Considering this, the lowest value of sulphur balance was determined in corn in treatment 
VI, with the difference between input and output of sulphur: -35.5 kg ha-1 S-SO4, while the 
highest level was determined for soybean in treatment VH, where the difference between 
input and output of sulphur was 335.5 kg ha-1 S-SO4. These could be explained because of 
high water shortage in 2011, when precipitation was 44 % below average (Zgorelec et al., 
2015), and as a result, there was no lysimeter or drainage leaching. Crops with low crop 
density (corn) had the highest losses through drainage (68.1 kg ha-1— VII) and lysimeter 
(251.3 kg ha-1— VII) water in comparison to crops a with high crop density (winter wheat, 
soybean, oilseed rape). 

Sulphur balance was analysed in the research of sulphate leaching in an organic crop 
rotation on sandy soil in Denmark, from 1994 to 1998. Considerable sulphate leaching 
occurred in the organic crop rotation, and it was equivalent to 60 % of the total input. The 
conclusion was that immediate sulphur deficiency may not occur, but negative sulphur 
balance must be expected in the long period (Eriksen, Askegaard 2000). A net negative total 
sulphur balance of 60 kg ha’ for sorghum, 54 kg ha” for rice and 68 kg ha’ for sugarcane, 
especially from non-sulphur applied plots were observed (Singh, 2014) in an experiment on 
sulphur balance in soils influenced by sulphur carriers and crops in India. 


Conclusions 

In this research, sulphur addition with FG, amelioration with zeolite tuff and CaCO; 
and fertilization with different nitrogen levels addresses important issues related to plant 
nutrition in especially intensive agriculture on less quality soil and the protection of the 
environment from pollution. Global S balance in the examined agroecosystem is in deficit 
except in the FG treatment. The results show that S application through FG is an acceptable 
and necessary measure for preventing soil degradation (in soils where S balance is in deficit 
in combination with calcification if soil pH is low). In this study S application resulted in 7.0 
% increase of S soil amount, with losses of 6.7 % through lysimeter water, and 18.3 % 
through drainage water. Long-term FG application could be beneficial however it is necessary 
to determine the optimal dose of FG for each soil type, climatic conditions and crop, and to 
find sustainable amount that would increase sulphur content in soil and minimize S leaching 
because excessive FG application can lead to excessive accumulation of S in soil, which 
could risk surface water and groundwater loads. 
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Abstract 

High soil quality is important for agricultural activities but pollution from cement 
industries is of great threat to sustainable agriculture in Nigeria. The aim of this study is to assess 
the effect of heavy metals in cement dust on the soil enzymes activities. 

This study investigated the activities of some selected soil enzymes, the microbial 
communities’ population, the concentrations of some heavy metals, some selected soil 
macronutrients and the physico-chemical properties (pH, moisture, and temperature) of the soil 
surface (0-15 cm), within the vicinity of a cement factory in Nigeria. 

The pH, total organic nitrogen (TON), total organic carbon (TOC) and heavy metals 
decreased with increasing distance from the cement factory. The microbial population around the 
factory was (34.00 + 2.64 CFU g’) and (92.00 + 1.00 CFU g’) at 100 m and the control site 
respectively. All the soil heavy metals contents analyzed are less than the permissible limit 
except Cd. Moreover, the enzyme activities increased away from cement factory except alkaline 
phosphatase that decreases with increase in distance. The enzymes activities (except alkaline 
phosphatase) correlated negatively with the content of all heavy metals, the pH, the TON, the 
TOC and temperature but positively correlated with moisture, phosphorous and microbial 
populations. The inverse correlations between most of the heavy metals, the pH and soil enzymes 
activities could be a result of the pollution. This means that the cement production exhibits a 
significant effect on the enzyme activities. 


Key Words: anthropogenic impact, bioaccumulation and biomagnification, 
geomagnification, soil enzymes, xenobiotics. 
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Introduction 

The biophysical environment contains definite amount of elements and compounds in a 
given proportion by nature and is considered polluted when there is an imbalance in the natural 
composition of all these compounds. The toxic metals, from the xenobiotics, entering the 
ecosystem may lead to geoaccumulation, bioaccumulation and biomagnifications in organisms in 
the environment (Ogunkunle and Fatoba, 2014). 

The xenobiotics can be released intentionally into the soil or dispersed by wind and rain 
to a far distance from the source of pollution. The concentrations of xenobiotics in soils are 
associated with biological and geochemical cycles and are influenced by anthropogenic 
activities, such as transport, waste disposal, industrialization, social and agricultural activities. 
They have an effect on environmental pollution and the global ecosystem. 

One industry that emits pollutants in the form of dust and gases which find their way into 
the soil is the cement industry (Addo, 2013). Dust from cement and other factories leads to 
considerable change in pH and accumulation of metals in the soil which may affect both the 
composition and physiological processes of microorganisms leading to a reduction in microbial 
biomass and enzyme activities (Utobo and Tewari, 2014). 

Soil enzymes are a group of enzymes that inhabits the soil and are continuously playing 
an important role in maintaining the soil biological, physical and chemical properties, fertility, 
and soil health (Das and Varma, 2011). All soils contain a group of enzymes that determine the 
soil metabolic processes (Baldrian, 2009) which in turn, depend on its physical, chemical, 
microbiological, and biochemical properties. 

Some of the soil enzymes including amylase, b-glucosidase, cellulose, chitinase, 
dehydrogenase, phosphatase, protease, and urease are released from plants and animals, organic 
compounds, and microorganisms in the soil. Soil enzymes are the mediators and catalysts of 
important soil functions and have a crucial role in carbon (f-glucosidase and f-galactosidase), 
nitrogen (urease), phosphorus (phosphatase), and sulphur (sulphatase) cycle (Karaca et al., 
2011). 

The importance of soil as an ecosystem cannot be underestimated. It serves as reservoir 
to a lot of nutrients and minerals which are essential for plant growth. However, soils around 
Lafarge cement factory in Sagamu are under threat of serious degradation as a result of cement 
dust deposition. It is therefore important to access the level of degradation for the purpose of 
conservation of biodiversity and protection of ecosystem around the factory. 

Also, to the best of our understanding, no work till date has reported the use of soil 
enzymes around Lafarge cement factory. Therefore, it is very important to assess the level of soil 
pollution around the study area using the enzyme activities, because the knowledge can be used 
as bioindicator of ecosystem perturbation, agricultural practices and xenobiotic pollution. 


Bulgarian Journal of Soil Science® 2018 Volume 3. Issue 2. www.bsss.bg 120 


Oludoye et al., 2018 





The aim of this study is to assess the effect of heavy metals in cement dust on soil 
enzymes activities around Lafarge cement factory in Sagamu, Sagamu Local Government Area, 
Ogun state, Nigeria. 

Materials and Methods 

The Study Area 

The study was carried out in the Lafarge cement factory which is situated in Sagamu, 
southwest Nigeria (6°50' 7°00' N; 3°45' 4°00' E). The factory is located 72 km southeast of 
Ibadan (largest city in West Africa), 32 km west of Ijebu Ode, 63 km southeast of Abeokuta and 
67 km northwest of Lagos. The factory was established, commissioned and became fully 
operational in 1978 with the present production capacity of 900 000 tonnes/year (Lafarge 
Cement WAPCO PLC, 2011). The area is characterized by high annual temperature, high 
rainfall, high evapotranspiration and high relative humidity which make it to be classified as 
humid tropical region (Ogunkunle and Fatoba, 2014). The soil type of Shagamu is Ferralitic and 
Ferruginous (Ilalokhoin et al., 2013). The climate is classified as humid tropical climatic zone 
and controlled by the Tropical Maritime and Tropical Continental air masses (Ilalokhoin et al., 
2013). The mean annual rainfall of Shagamu for 2015 was 1 100 mm (NIMET, 2015). 





Figure 1. Google map of Lafarge factory 


Collection of Samples 

Sampling was carried out in the wet seasons of 2015 in the area surrounding the factory 
excluding the northern part of the factory because of the presence of swamps. At each sampling 
location, three subsamples were collected at the depth of 0 — 15 cm within 10 m radius of a point 
using a soil auger. The samples were bulked together to obtain a composite sample and recorded 
against the coordinate position. A total of 15 composite topsoil samples were collected in the 
west, the south and the east parts at 100 m, 200 m, 300 m, 500 m and 1000 m distance from the 
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factory respectively. The coordinates at every sampling points were taken by a hand-held GPS 
(Garmin 2H). 

Control samples were collected at a distance of about 15 km from the cement factory 
where there was no record of cement dust pollution. The soil samples collected were placed in 
polyethylene bags, properly labeled and transported to the laboratory. Before chemical analysis 
in the laboratory, the soil samples were air—dried, sieved through a 2—mm sieve mesh to remove 
coarse material and debris, and pulverized into powder. The soil was transferred to the storage 
room and was stored at 30 °C until the time of analysis. Microbial and enzyme analysis were 
done within 24 to 48 h. 


Soil Physico-Chemical Analysis 

The soil pH was measured in a soil water suspension (1 : 2, soil : water) using glass 
electrode and temperature was also measured by soil thermometer (Frederick et al., 2014) while 
the moisture content was determined gravimetrically after drying at 105 °C according to AOAC 
(2005). The following soil nutrients were selected and analyzed as follows: total organic carbon 
was determined based on modified Walkley-Black chromic acid wet oxidation method according 
to Sato et al. (2014). Total phosphorus (TP) was determined according to Adesanwo et al. 
(2013). Total nitrogen was determined by Kjeldahl method according to AOAC (2005). 


Heavy Metal Analysis in Soil 
The soil samples (air-dried and digested) were analyzed for heavy metals content using Peking 
Elmer Model 403 atomic absorption spectrophotometer (AAS) according to AOAC (2005) after 
nitric acid digestion of the soil. The heavy metals studied were copper (Cu), cobalt (Co), lead 
(Pb), manganese (Mn), cadmium (Cd), nickel (Ni) and zinc (Zn). In addition, calctum (Ca) and 
aluminum (Al) content of the soil samples were also analyzed because they are one of the main 
constituent elements in the cement. 


Enzymatic Analysis and Microbial Population in Soil 

Enzyme substrates were purchased from Sigma-Aldrich Inc., England. Enzyme activities 
were measured according increase of resultant or decrease of substrate by colorimetric 
determination methods. Soil urease, acid phosphatase, alkaline phosphatase and dehydrogenase 
were determined according to Ofoegbu ef al. (2013). B-glucosidase activity was determined 
according to Ramamoorthy (2015) and Zhang (2010). Soil invertase activity was determined 
according to the method described by Frankenberger Jr. and Johanson (1983) cited by Effron et 
al. (2004). Cellulase activity was measured with a method adapted from the one described by 
Alef and Nannipieri (1995). The method employed for determining soil amylase activity was 
essentially based on starch hydrolysis according to Khan et al. (2007). 

Also, the number of microbial population was determined by the method of serial 
dilutions (Kulandaivel et al., 2015). Soil samples were subjected to serial dilution and spread on 
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nutrient agar plates and potato agar plates, incubated for 48 h to grow the microbial colonies 
proper. Colony forming units were counted by colony counter. 

Statistical Analysis 

Data were subjected to two-way Analysis of Variance (ANOVA) and Dunnett’s multiple 
comparison tests at p < 0.05 probability level using Statistical Package for the Social Sciences 
(SPSS) version 20. Spearman correlation analysis was also employed to determine the 
relationships existing among the heavy metals content, the pH, the microbial population, enzyme 
activities and soil nutrient contents. 


Results and Discussion 

The results from the Analysis of Variance (ANOVA) showed that there were significance 
mean differences at (p < 0.05) between the soil samples taken at various distances and the con- 
trol sample for all the enzyme activities, heavy metals; the soil nutrients except total organic 
nitrogen (f = 0.18), the soil physiochemical properties except pH values and lastly significant for 
all the microbial population constituents.Fang et al. (2017) analyzed the activity of a number 
of enzymes, in three types of land (farmland, woodland, and grassland) to evaluate soil pollution 
by heavy metals (Pb, Zn, and Cd). Their results showed that urease was the most sensitive 
enzyme to heavy metal pollution in the farmland. This is in agreement with our results where 
urease had the lowest enzyme activities (at 100 m). Caravaca et al. (2017) also reported that 
Urease activity was only decreased in the zone nearest to cement factory; and glucosidase 
activities declined significantly with decreasing distance to the cement factory in comparison 
with that in the unaffected zone which is in agreement with our work. Furthermore, Chen et al. 
(2017) reported that the soil moisture, contents of silt and total nitrogen were the most important 
variables responsible for soil extracellular enzyme activity variation. 


Enzyme activities and the physicochemical properties of the soil 

There was significant difference (p < 0.05) in soil enzymatic activities with respect to the 
distance from the factory except alkaline phosphatase. This implies that the enzymatic inhibition 
could originate from the factory (cement dust). The highest activity was recorded for invertase 
(517.33 + 2.08 ug g-1) at the control site while the lowest activity was recorded for urease (0.14 
+ 0.03 pg g-1) at 100 m distance respectively. The enzyme inhibition could be attributed to the 
binding of the heavy metals at the active site of the enzymes which could subsequently inhibit 
their catalysis. According to Poli et al. (2009) heavy metals is a strong enzyme inhibitors. 

The pH range between 6.62 + 0.10 (at the control) to 10.77 + 0.81 (at 100 m) indicated 
that the factory soil was alkaline. The alkalinity was due to deposition and liming by calcium ion 
(Ca2+) from cement dust. Invertase had the highest enzymatic activity (517.33 + 2.08 ug g-1) at 
pH of 6.62 + 0.10 (at the control) while the lowest value was recorded for urease (0.14 + 0.03 pg 
g-1) at pH of 10.77 + 0.81. This may be as a result of higher sensitivity of urease to alkaline pH 
as similarly reported by Frederick et al. (2014). Generally, the lower enzymes activities could be 
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due high pH around the factory as most of these enzymes could only function optimally at lower 


pH except alkaline phosphatase. 


Table 1. The relationship between the enzyme activity and the sampling distance from different 






































locations 

Enzyme Distance - radius (m) 

(ug g-1) | 100 200 300 500 1000 15000 

URS 0.14 +0.03 | 0.21+0.04 0.30+0.06 0.37+0.06 0.48+0.08 0.97+0.01 

ALP 68.67+4.16 | 52.6745.03 | 57.67+4.04 54.33+42.08 32.6741.57 | 26.67+0.58 

ACP 101.3343.0 | 114.67+4.16 | 134.3345.51 | 145.00+4.00 | 163.00+5.0 | 177.67+1.5 
6 0 3 

BGLU _ | 2.47+0.38 | 5.17+0.25 6.74+0.31 7.83+0.25 9.27+0.32 11.23+0.25 

DHA 1.60+0.3 2.27+0.60 3.0+0.79 3.63+0.75 4.541.10 7.7740.15 

INVS 120.67422. | 156.33420.4 | 212.33436.6 | 268.33469.4 | 349466.26 | 517.3342.0 
75 0 7 0 8 

CLS 17.6743.06 | 23.00+4.58 | 28.67+4.50 33.3343.51 40.6743.51 | 6542.64 

AMY 4.7140.95 | 6.21+1.73 7.49+1.48 8.76+1.19 10.7840.58 | 12.17+0.12 























Values are mean + standard error of triplicate values. 

KEY: Urease (URS) (ug N-urea gary soi 5h’'), Alkaline phosphatase (ALP) (ugp-nitrofenol g” ary 
sot h’'), Acid phosphatase (ACP) (ugp-nitrofenol gary sou h'), B-glucosidase (BGLU) (ugp- 
nitrofenol go ary soil h”), Dehydrogenase (DHA) (ugtriphenylformazan gary soil 24h"), Invertase 
(INVS) (ug reducing sugar So ary soil 24h”), Cellulase (CLS) (ug glucose @ ary soil 24h”), Amylase 
(AMY)(ug starch degraded —_o soit 2-5h""), 


The study showed that the soil pH around factory was positively correlated for alkaline 
phosphatase while it was negatively correlated with all other enzymes, because the enzyme 
requires alkaline pH for optimal activity. This is in agreement with the work of Karaca et al. 
(2011) and Caravaca et al. (2017). 

Also, there were positive correlations between the soil moisture content and the enzymes 
activities. This showed that moisture is required for optimum enzyme activity (through high 
microbial activities) in the soil. Also, high moisture has been predicted to support high microbial 
activities in the soil (Ahmed et al., 2018). 

There was significance difference (p < 0.05) between the mean values of the pH, soil 
moisture (except the temperature) and the enzyme activities at different sampling sites compared 
with the control samples. 


The heavy metals, Calcium, Aluminum and the enzyme activities 
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The heavy metal concentrations decrease with increasing distance from the factory. The 
lowest concentration was observed at the control soil sample (15000 m) except zinc. The 
concentration of zinc at the control site (15.32 mg kg’) was higher than the concentrations of Zn 
around the factory. There were statistically significant variations in the concentrations of Cu, Pb, 
Mn, Cd, Ni, Co, Zn, Al and Ca at some distances, when compared to the control samples (p < 
0.05). This implied that heavy metals are from the cement dust pollution. Moreover, Ca has the 
highest concentration at all sampling points which are also all statistically significant to the 
control and could be the main contributor to soil pollution in the environment. This is because Ca 
is an important element in the cement kiln process. 

Dunnett’s Multiple Comparison tests also revealed further that for Cu and Mn, there were 
no significant differences in the mean soil samples at 500 m, 1000 m compared to the control. 
Cadmium showed significant differences at different locations except at 1000 m only. We 
reported previously that Zn showed no significant difference at 200 m, 300 m, 500 m and 1000 
m compared to the control and Al at 1000 m to the control samples (Oludoye and Ogunyebi, 
2017). 

There were also negative correlations between the enzyme activities and heavy metal 
concentrations around the factory except the alkaline phosphatase that was positively correlated. 
Heavy metals have been reported to inhibit enzymes singly or in combination by synergism 
(Frederick et al., 2014). Therefore, this implied that some of these metals might inhibit the 
activities of all the enzymes studied except alkaline phosphatase activities which are also in 
accordance with the work of Das and Varma (2011). 

In order to determine the magnitude of the pollution, the heavy metals content was 
compared with Nigeria Federal Environmental Protection Agency (FEPA) maximum permissible 
concentration (MPC) for soil. It was observed that the concentrations of all the heavy metals 
except cadmium are within the MPC of FEPA. The soil Cd concentration is higher than the MPC 
at all distances except at the control. This shows that Cd concentration is very high around the 
cement factory. Cadmium is one of the trace elements found in raw mill and coal which are used 
in the formation of clinker (the most important constituent of cement) (Cipurkovic et al., 2014). 
This high Cd concentration is unhealthy for the environment and human health because there are 
some farming activities currently going on in the vicinity of the cement factory. Cadmium could 
subsequently enter the food chains through the crops from the farms. 

Table 2. FEPA and WHO maximum permissible concentrations for soil 























Heavy metals WHO limit (mg/kg) FEPA limit (mg/kg) 
Zn 50 300-400 

Cu 36 70-80 

Pb 85 1.6 

Cd 0.8 3.0 

Mn Not fixed Not fixed 
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Ni 35 18 








Co Not fixed 20 











However, it is not only cement production that could release Cd into the environment but 
geologic processes (ferratillisation) could also contribute to the high concentration around the 
factory. It was also reported that high soil pH caused by limestone deposition (a component of 
cement) could immobilize Cd in the soil (Yun and Yu, 2015). Furthermore, Kim (2014) reported 
that Cd is maximally adsorbed at high pH and high temperature which significantly affected the 
adsorption capacity in alkaline soil. 

Table 3. The mean value of heavy metal concentrations in the soil around the factory 






































Heavy Average distance - radius (m) from east, west and south part of the factory 
metal 100 200 300 500 1000 15000 

(mg kg”) 

Cu 6.52 + 0.80 | 5.69 + 0.50 4.62 + 0.79 3.71 £0.58 | 3.04+0.53 2.16 + 0.02 
Pb 1.3440.06 | 1.24+0.07 1.14+ 0.09 1.06 +0.08 | 0.96+0.09 0.27 + 0.01 
Cd 9.27+0.42 | 8.63+0.55 7.58 + 0.44 6.60 0.59 | 5.7+0.42 2.56 + 0.09 
Mn 9.87+0.19 | 9.25+0.24 8.42 + 0.49 8.00 0.36 | 7.38+0.27 6.77 + 0.04 
Ni 3.52 + 0.09 | 3.30+0.25 3.16 +0.21 2.97 0.17 | 2.86+0.16 2.09 + 0.07 
Co 16.97 + 0.72 | 16.2240.61 | 15.05+0.42 | 14.14 +0.54 | 13.7440.81 | 10.48 + 0.10 
Zn 17.25 + 0.21 | 16.68+0.26 | 16.15+40.37 | 15.47 40.70 | 14.7140.54 | 15.32 + 0.25 




















Values are mean + standard error of triplicate values 


Table 4. The mean value of Ca and Al concentrations in the soil around the factory 

















Metal Average distance - radius (m) from east, west and south part of the factory 

(mg kg") | 100 200 300 500 1000 15000 

Ca 201.05412.0 | 160.924 144.02411.5 | 133.80+9.8 | 97.0641.53 | 84.13 + 0.26 
11.08 8 6 

Al 4.15+0.06 | 3.65+0.34 3.17+40.32 2.5340.45 = | 2.1340.17 1.81 + 0.03 























The soil nutrients analysis 
There was decrease in the concentration of available phosphorus at various distances 








which showed a significant mean difference to the control at 15000 m. This can attributed to the 
precipitation of phosphorus into calcium phosphate at high pH (Cao et al., 2007) while there was 
gradual decrease in the total organic carbon and nitrogen content of the soil from the factory to 
the control site, respectively. 

Dunnett’s Multiple Comparison tests showed that for there were no significant differences be- 
tween the means for all the soil samples taken at varying distances and the control sample, for 
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organic nitrogen. In addition, there was no significant mean difference at 300 m, 500 m and 1000 


m for total organic carbon. 


Table 5. The mean values of selected soil nutrients with respect to distance 






































Distance (m) Phosphorous (part per_ | Nitrogen (mg kg”) Carbon (%) 
million) 
100 21.43 + 2.59 178.00 + 4.00 14.27 + 1.02 
200 24.57 + 2.61 164.00 + 4.00 12.53 + 2.06 
300 27.60 + 3.53 144.00 + 3.00 9.44 + 0.70 
500 29.47 + 3.19 133.00 + 3.00 7.59 + 0.78 
1000 31.03 + 2.06 110.00 + 2.00 6.72 + 0.44 
Control 58.30 + 0.18 98.00 + 4.00 6.30 + 0.20 
Values are mean + standard error of triplicate values. 
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Figure 2. The relationship between the pH, moisture, and temperature with respect to distance 


Soil microbial population 
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The microbial communities of the soil studied were significantly affected as a result 
of cement dust deposition. The result showed that the microbial communities of the control were 
statistically significantly higher (p < 0.05) than the polluted soil around the factory, which is in 
agreement with the work of Kulandaivel et al. (2015). Also, Caracava et al. (2017) reported that 
the prolonged exposure of soil to cement dust could result in shifts in the soil microbial function. 
This could also be attributed to the decrease in microbial biomass due to reduced enzymes 
activity.Spearman correlation analysis revealed that there was positive correlation between the 
microbial communities and the enzyme activities. This implied that enzymes are primarily 
derived from microorganisms but can also originate from plants and animals. Hence, factors 
affecting microbial activities will invariably affect soil enzymes. Moreover, the microbial count 
is also negatively correlated with the soil pH. This could be due to the fact that some of the 
microbes could not survive in alkaline conditions as this could affect the osmotic balance of their 
cells and consequently, cause reduction in their cell biomass. 
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Figure 3. The relationship of microbial population from different soil samples with respect to 
distance from the factory. 
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Table 6. The Spearman's correlation coefficient among the enzyme activities and the heavy 


metal concentrations 


Pe a NT co ee | 





* Correlation is significant at the 0.05 level (2-tailed). 
** Correlation is significant at the 0.01 level (2-tailed). 


Conclusion 

All the heavy metals content of soil around the cement factory are below the MPC of 
FEPA except cadmium. The high Cd content could have emanated from the factory cement dust 
because cadmium is one of the trace elements in the fuel that is used for kiln production. Also, 
soil geologic processes and high pH around the factory could also result in Cd immobilization in 
the soil. Moreover, there is low nutrients content (phosphorus), low microbial population and 
high pH within the vicinity of the cement factory. All these factors had a significant adverse 
effect on and are inversely correlated with the soil enzymes. This indicates the ability of soil 
enzymes early soil pollution signal (or indicator) despite the permissible level of heavy metals 
(except cadmium). 
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The heavy metals content (especially cadmium) in the top soil (0-15 cm) is within 
the reach of plant roots and could have negative effects on the natives who are already using the 
sites around the factory as agricultural farms for planting food crops such as sugar cane, cassava, 
banana and maize respectively. 

Therefore, there should be comprehensive ecological risk assessment of the soil 
ecosystem around the factory. In addition, human health risk assessment should be carried out to 
determine carcinogenic and non-carcinogenic risk through dietary exposure to some of these 
heavy metals (especially cadmium). This would help prevent the pollutants (heavy metals) from 
entering the food chain and help reduce the exposure of the inhabitants and the farmers. Also, 
farming activities should be totally discouraged around the factory. 


Reference 

Addo, M.A., Darko, E.O., Gordon, C., Nyarko, B.J.B., Gbadago, J.K., Nyarko, E., 
Affum, H.A., Botwe, B.O., 2012. Evaluation of heavy metals contamination of soil and 
vegetation in the vicinity of a cement factory in the Volta region, Ghana. International Journal 
of Science and Technology 2(1): 40—S0. 

Adesanwo, O. O., Ige, D. V., Thibault, L., Flaten D., Akinremi, W. 2013. Comparison of 
Colorimetric and ICP Methods of Phosphorus Determination in Soil Extracts. Communications 
in Soil Science and Plant Analysis 44(21): 3061-3075. Doi: 
https://doi.org/10.1080/00103624.2013.832771 . 

Ahmed, M. A., Sanaullah, M., Blagodatskaya, E., Mason-Jones, K., Jawad, H., 
Kuzyakov, Y., Dippold, M. A. 2018. Soil microorganisms exhibit enzymatic and priming 
response to root mucilage under drought. Soil Biology & Biochemistry 116: 410-418. Doi: 
https://doi.org/ 10.1016/j.soilbio.2017.10.041 . 

Alef, K., Nannipieri, P. 1995. Cellulase activity. In: Alef P (ed) Methods in applied soil 
microbiology and biochemistry. Academic Press, London, pp 382-384. 

AOAC. 2005. Official methods of analysis. Association of official analytical chemists: 
Agricultural Chemicals; Contaminants; Drugs. 18" Ed. Association of Official Analytical 
Chemists Inc. Arlington, Virginia, USA. 

Baldrian, P., 2009. Microbial enzyme-catalyzed processes in soils and their analysis. 
Plant and Soil Environment 55 (9): 370-378. 

Bradford, M.M., 1976. A rapid and sensitive method for the quantitation of microgram 
quantities of protein utilizing the principle of protein dye binding. Analytical Biochemistry 72(1- 
2): 248-254. Doi: https://doi.org/10.1016/0003-2697(76)90527-3 . 

Cao,X., Harris, W. G., Josan, M. S. and Nair,V. D., 2007. Inhibition of calcium 
phosphate precipitation under environmentally-relevant conditions. Science of the Total 
Environment. 383(1-3): 205-215. 





Bulgarian Journal of Soil Science® 2018 Volume 3. Issue 2. www.bsss.bg 130 


Oludoye et al., 2018 





Chen, H., Luo, P., Wen, L., Yang, L. Q., Wang, K. L., Li, D. J. 2017. Determinants of 
soil extracellular enzyme activity in a karst region, southwest China. European Journal of Soil 
Biology 80: 69-76. Doi: https://doi.org/10.1016/j.ejsobi.2017.05.001 . 

Das, S.K., Varma, A., 2011. Role of enzymes in maintaining soil health. In: In: Soil 
Enzymology, Soil Biology 22. Shukla, G., Varma, A. (Eds.). Springer-Verlag, Berlin Heidelberg, 
pp.25-42. 

Cipurkovic, A., Trumic, I., Hodzic, Z. SelimbaSic, V. and Djozic, A. (2014). Distribution 
of heavy metals in Portland cement production process. 5(6): 252-259. 

Effron, D., de la Horra, A.M., Defrieri, R.L., Fontanive, V., Palma, P.M., 2004. Effect of 
cadmium, copper, and lead on different enzyme activities in a native forest soil. Communication 
of Soil Science and Plant Analysis 35(9-10): 1309-1321. Doi: https://doi.org/10.1081/CSS- 
120037548 . 

Egejuru, L.O., Akubugwo, E.I., Ude, V.C., Ugbogu, O.C., Ugbogu, A.E., 2014. 
Evaluation of physicochemical properties, microbial loads and enzymes activity studies of 
agrochemicals on the Imo-River basin farm in Imo State Nigeria. International Journal of 
Current Microbiology and Applied Sciences 3(2): 776-786. 

Fang, L. C., Liu, Y. Q., Tian, H. X., Chen, H. S., Wang, Y. Q., Huang, M. 2017. Proper 
land use for heavy metal-polluted soil based on enzyme activity analysis around a Pb-Zn mine in 
Feng County, China. Environmental Science and Pollution Research 24(36): 28152-28164.Doi: 
https://doi.org/ 10.1007/s11356-017-0308-4 . 

Frankenberger Jr., W. T., Johanson, J. B., 1983. Method of measuring invertase activity 
in soils. Plant and Soil 74: 301-311. 

Frederick, O.0., Iroha, A.E., Oswald, E.C., 2014. Evaluation of the concentration of 
selected heavy metals and the effects on soil enzymatic activities in an Abandoned Cement 





Factory in NigercemNkalagu and its environs. International Journal of Biochemistry Research 
and Review 4(1): 16-27. Doi: https://doi.org/10.9734/IJBCRR/2014/5019 . 

Tlalokhoin, P.O., Otaru, A.J., Odigure, J.O., Abdulkareem, A.S., Okafor, J.O., 2013. 
Environmental impact assessment of a proposed cement plant in southwestern Nigeria. 
IOSRJournal of Environmental Science, Toxicology and Food Technology 3(5): 83-99. 

Karaca, A., Cema, C.C., Turgay, O.C., Kizilkaya, R., 2011. Soil enzymes as indicator of 
soil quality. In: Soil Enzymology, Soil Biology 22. Shukla, G., Varma, A. (Eds.). Springer- 
Verlag, Berlin Heidelberg, pp.119-148. 

Khan, S., Cao, Q., Hesham, A.E., Xia, Y., He, J., 2007. Soil enzymatic activities and 
microbial community structure with different application rates of Cd and Pb. Journal of 
Environmental Science 19(7): 834-840. Doi: https://doi.org/10.1016/S1001-0742(07)60139-9 

Kim, M. (2014). A study on the adsorption characteristics of cadmium and zinc onto 
acidic and alkaline soils. Environmental Earth Sciences. 72, (10): 3981-3990. 


Bulgarian Journal of Soil Science® 2018 Volume 3. Issue 2. www.bsss.bg 131 


Oludoye et al., 2018 





Kulandaivel, S., Nagarajan, S., Priyanga, A., Saravanapandian, R., Thangarani, A., 2015. 
Effect of cement dust pollution on microbial properties and alkaline phosphatase enzyme activity 
in soil. International Journal of Current Microbiology and Applied Science4(2): 641-649. 

Lafarge Cement WAPCO PLC. 2011. Lafarge WAPCO: Supporting Nigeria’s 
development, Annual Report—Lafarge Cement WAPCO PLC. pp. 2-4. 

Nayak, D.C., Sen, T.K., Chamuah, G.S., Sehgal, J.L., 1996. Nature of soil acidity in 
some soils of Mainpur. Journal of Indian Society of Soil Science 44(2): 209-212. 

Nelson, N. 1944. A photometric adaptation of the Somogyi method for determination of 
glucose. Journal of Biological Chemistry 153: 375-380. 

NIMET. 2015. Seasonal rainfall prediction and socio-economic implications for Nigeria. 
Federal Government of Nigeria, Nigerian Meteorological Agency (NIMET), Abuja, Nigeria, 
Available at [access date : 30.11.2017]: http://nimet.gov.ng/sites/default/files/publications/2012- 
seasonal-rainfall-prediction.pdf. 

Ofoegbu, C.J., Akubugwo, E.I., Dike, C.C., Maduka, H.C.C., Ugwu, C.E., Obasi, N.A., 
2013. Effects of heavy metals on soil enzymatic activities in the Ishiagu mining area of Ebonyi 
State-Nigeria. JOSR Journal of Environmental Science, Toxicology and Food Science 5(6): 66- 
71. 








Ogunkunle, C.O., Fatoba, P.O., 2014. Contamination and spatial distribution of heavy 
metals in topsoil surrounding a mega cement factory. Atmospheric Pollution Research 5(2): 270— 
282. Doi: https://doi.org/10.5094/APR.2014.033 . 

Oludoye, O. O. and Ogunyebi, L. A., 2017. Nutrients Assessment of Tropical Soils 
Around A MegaCement Factory In Southwest Nigeria. Journal of Ecological Engineering. 18(2): 
18-25. 

Poli, A., Salerno, A., Laezza, G., di Donato, P. and Dumontet, S., 2009. Heavy metal 
resistance of some thermophiles: potential use of alpha-amylase 





from Anoxybacillusamylolyticus as a microbial enzymatic bioassay. Research in Microbiology. 
160: 99-106. 

Ramamoorthy, D., 2015. Assessment of heavy metal pollution and its impacts on soil 
physical, chemical properties and f-glucosidase activities in agricultural lands, Puducherry 
region. Project Report. Department of Science, Technology and Environment Government of 
Puducherry. 45pp 

Sato, J.H., de Figueiredo, C.C., Marchao, R.L., Madari, B.E., Benedito, L.E.C., Busato, 
J.G., de Souza, D.M. 2014. Methods of soil organic carbon determination in Brazilian savannah 
soils. Scientia Agricola 71(4): 302-308. Doi: http://dx.doi.org/10.1590/0103-9016-2013-0306 . 

Utobo, E.B., Tewari, L., 2014. Soil enzymes as bioindicators of soil ecosystem status. 
Applied Ecology and Environmental Research 13(1): 147-169. Doi: 
https://doi.org/10.15666/aeer/1301_ 147169 . 











Bulgarian Journal of Soil Science® 2018 Volume 3. Issue 2. www.bsss.bg 132 


Oludoye et al., 2018 





Wainwright, M., Duddridge, J. E., Killham, K., 1981. Assay of a-amylase in soil and 
river sediments: its use to determine the effects of heavy metals on starch degradation. Enzyme 
and Microbial Technology 4: 32-34. 

Wang, Y.P., Shi, J.Y., Wang, H., Lin, Q., Chen, X.C., Chen, Y.X., 2007. The influence of 
soil heavy metals pollution on soil microbial biomass, enzyme activity, and community 
composition near a copper smelter. Ecotoxicology and Environmental Safety 67(1): 75—81. Doi: 
https://doi.org/10.1016/j.ecoenv.2006.03.007 . 

Wright, A.L., Reddy, K.R., 2001. Heterotrophic microbial activity in northern Everglades 
wetland soils. Soil Science Society of America Journal 65(6): 1856-1864. Doi: 
https://doi.org/10.2136/sssaj2001.1856. 

Yun, S. and Yu, C., 2015. Immobilization of Cd, Zn, and Pb from Soil Treated by 
Limestone with Variation of pH Using a Column Test. Journal of Chemistry. 
http://dx.doi.org/10.1155/2015/641415 . 

Zhang, X., Li, F., Liu, T., Xu, C., Duan, D., Peng, C., Zhu, S., Shi, J., 2013. The 
variations in the soil enzyme activity, protein expression, microbial biomass, and community 
structure of soil contaminated by heavy metals. ISRN Soil Science Article ID 803150. Doi: 
http://dx.doi.org/10.1155/2013/803150 . 








Bulgarian Journal of Soil Science® 2018 Volume 3. Issue 2. www.bsss.bg 133 





Opportunities for growing lavender in areas facing 
natural constraints exemplified by the 
Borika village, Bulgaria. 

Zorintsa Mitreva, Veselin Pankov 





N. Poushkarov Institute of Soil Science, Agrotechnologies and Plant Protection 
1331 Sofia, 7 Shosse Bankya str. 
Corresponding Author: Zorintsa Mitreva, e-mail: zuza_na@aby.bg 





Abstract: 

Areas facing natural constraints represent a serious obstacle for Bulgarian farmers and 
to a great extent limit their development. In this regard, it is necessary to look for alternative 
options for optimal use of these areas so that they are as efficient as possible despite the 
limitations they have. 

The paper presents the basic requirements of lavender in terms of climatic and soil 
conditions. There is also a soil and climate characterization of a village located in area of high 
altitude - the land of the village of Borika, Sofia in order to be able to indicate which soil 
types in the Borika land are best suited for lavender growing. 


Key words: field ratings, soil, soil characteristics, lavender 


Introduction 

The land of the village of Borika is located in the Southwest region of Bulgaria, 
Thtiman municipality, Sofia district. It falls in a mountainous area with an altitude of 700 - 
990 m. The total area of Borika is 294,75 ha. 

The soil resources in the area were studied in the large-scale soil mapping with scale 
1:10 000, conducted by scientists of Nikola Pushkarov Institute of soil science, 
agrotechnologies and plant protection. In the land of the village there are 11 soil types, for 
each are presented detailed soil characteristics, which are related to the cultivation of 
lavender: soil texture (contains soil texture of the plowed horizon with < 0,01 mm particle 
size in %, the soil texture of the soil under the plowed horizon < 0,01 mm particle size in % 
(Kachinskiy, 1958), humus horizon depth (cm), soil profile depth (cm), soil reaction in HO 
(pH in water extraction — measured potentiometrically (WTW 720 pH meter)), texture 
coefficient (which is the % clay content in the B horizon divided by the % clay content in the 
A horizon with particles < 0.001 mm (Kachinskiy, 1958), humus content was estimated by the 
Turin method (Donov et al., 1974). 

According to the European Union Regulation currently (2018) in force (Regulation 
(EU) No 1305/2013) Annex III the biophysical criteria for the delimitation of areas facing 
natural constraints Borika village falls into an area facing natural constraints. 
Borika village is a mountainous region that has soil types with unfavorable texture (low 
content of particles <0,01 mm) and stoniness, shallow rooting depth (depth<30 cm) and poor 
chemical properties (acidity<S5). 
Lavender 

Lavender (Lavandula vera L), belongs to the Lamiaceae family, which is major oil 
crop for the temperate climate zone and for Bulgaria. The Mediterranean countries, France, 
Spain, Italy and the Balkans are where the lavender originated (Reed, 2000). Lavender is a 
good honey-bearing and soil-strengthening plant. It is a perennial shrub, evergreen tuft. 
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Lavendula angustifolia grows 60-80 cm and has fragrant violet flower spikes of 3-7 
cm (Reed, 2000). Spike lavender or L. latifolia (syn. L. spica) is also native to the western 
Mediterranean region and Portugal. It is later flowering and shorter than Lavendula 
angustifolia and the oil has a distinct camphoraceous odor (Reed, 2000). Important species are 
Lavendula dentata and Lavendula stoechas (Ernst, 2017). The root reaches a depth of 2-3 m 
(Koedzhikov et. al., 1971). It requires warm weather, strong solar radiation and lack of strong 
air currents. These requirements are largely satisfied on southern mountain slopes, protected 
from winds and over 800-1000 m height above the sea level (Koedzhikov et. al., 1971). 
During flowering changing and cloudy weather, as well as frequent precipitation, have the 
effect of lowering essential oil yields and deteriorating quality (Koedzhikov et. al., 1971). 

Lavender tolerates low temperatures at an early stage of development of - 8°, - 10°C 
(Koedzhikov et. al., 1971). Adult plants with sufficient snow cover can withstand up to - 25°, 
- 30°C (Koedzhikov et. al., 1971). High summer temperatures adversely affect oil quality 
(Reed, 2000). 

Lavender is very demanding in terms of light. When there is strong shading the oil in 
the flower decreases by 50% (Koedzhikov et. al., 1971). Requirement for moisture is not so 
demanding, it is drought-resistant, and excessive moisture gives negative results (Koedzhikov 
et. al., 1971). High summer temperatures adversely affect oil quality (Reid, 2000). 

The crop is not demanding with respect to soils. It develops optimally on loose, air- 
permeable, easily absorbing heat and rich in humus and lime soils (Koedzhikov et. al., 1971). 
Of particular importance to the high productivity of lavender is the soil reaction. Soils of 
neutral (7) pH are preferred but successful crops have been produced over the range 5.8 to 8.3 
pH (Reed, 2000). Lavender is best established on sandy loamy soils (Ernst, 2017). The crop is 
most successful when grown on well-drained soils (Ernst, 2017). Inappropriate are heavy, 
waterlogged and acidic soils with a high level of groundwater, as well as soils with high clay 
content (Koedzhikov et. al., 1971). 

The main objectives of the study are: 

(1) Present soil and climatic characterization of the land. 

(2) Consider the lavender requirements for the main soil climatic characteristics to 
indicate which soil types in the land are suited for cultivation and which are impossible for 
cultivation. 


Materials and Methods 
The data, interpreted below are from the Vakarel and Sofia meteorological stations, 
published (officially) by NIMH-BAS and summarized for a minimum of 45 years (Climate 
reference book for Bulgaria, 1978, 1983, 1990).Of particular importance in agriculture is the 
duration of the period with average daily air temperatures above a certain limit and the 
temperature sums for it. "Methods for working on the agricultural land cadastre in the 
Republic of Bulgaria" (Petrov et al., 1988) uses these characteristics with a lower limit of 10 
°C. The temperature sums for this period are called "active temperature sums" and are defined 
by Equation 1: 
YT °C> 10°C= (ty +t + ty +... + ty) (1) 
Where: 
> T °C> 10°C- the sum of effective temperatures (for the period with 
average daily temperature> 10°C) in °C; 
ti, to, tz, ... ty - consecutive observations of mean 24-hour temperatures 
of air in “C; 
1, 2, ... n - index for serial number of the day through 
the period of measurement. 
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The land capability approach uses the so-called limitation method, which implies 
comparing land characteristics of a land unit with critical limits of each capability class. The 
values of land characteristics are tested first against the limits of the highest capability class 
and, unless all the limits are met, the land unit automatically falls to the next lower class (Van 
Lanen, 1991). This sieving process continues until a capability class is found at which all the 
limits are met. This means that if one limitation is of sufficient severity to lower the land unit 
to a lower class, it is allocated to that class, no matter how favorable the other land 
characteristics may be (law of minimum) (Van Lanen, 1991). The results are harmonized with 
classes from the FAO framework (FAO, 1976, 2007). 

The nature of the severest limitation, e.g. erosion, excess water or climatic limitation, 
is specified in the capability assessment (Van Lanen, 1991). Refinements are possible by 
making the suitability class ratings dependent on more than one limiting land characteristic 
(Rosa, Diepen, 2002). 

Simple limitation method compares the plant requirements with its corresponding 
qualitative land and climatic characteristics and the most limiting characteristics defines land 
suitability class (Darvishi-Foshtomi et al., 2011). 

The limitation method concerns the selection of the most restricting criterion rate and 
considering it as the total score for a land unit (Sharififar, 2012). This method defines land 
classes with regard to number and intensity of limitations (Chadar et al., 2018, Sys et al., 
1991). 


Results and Discussion 


Climatic characteristic 

The land of the village of Borika falls into the Moderate Continental Climatic Sub- 
Area and more precisely in the climatic region of the high fields of Western Central Bulgaria. 
The climatic region of the highlands of Western Central Bulgaria consists of separate valleys 
located in the western part of Central Bulgaria and surrounded by many low and medium 
altitude mountains. Winter in the region is cold, but the summer is cooler. Typical of the area 
are late spring and early autumn frosts, due to the valley nature of the terrain. 


Precipitation 


Table 1. Average monthly and average annual precipitation. (Vakarel meteorological station) 





Month I | Wl} WV} V |} Vi} Vit | Vol | i | X | XI} XI | Annual sum 















































mm/m? | 42 | 37 | 44 | 55 | 85 | 82 | 63 44 40 | 49 | 49 | 47 635 





Table 1. shows that the maximum rainfall is in June and the minimum is in February 
in the nearest Vakarel meteorological station. The average annual precipitation is 635 mm / 
m’ (Climate reference book for Bulgaria, 1990). 

The mountainous region of Borika has sufficient precipitation for lavender as the crop is not 
very demanding in terms of water needs. 


Temperature characteristics 

As a start date of sustained air temperature above 10°C for the area is indicated 29 
April and for the final is 9 October (Climate reference book for Bulgaria, 1983, p.173). The 
duration of days with sustained air temperature > 10°C for the area is 162 days and the 
temperature sum is 2510°C for this period. (Climate reference book for Bulgaria, 1983 (p. 
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179)) as is indicated by the nearest Vakarel meteorological station. This implies good 
conditions for the development of lavender during vegetation. 


Characteristics of solar radiation 

One of the most important condition for lavender growing is the amount of solar 
radiation, the hourly and daily sums of direct photosynthetic active radiation in a cloudless 
sky on the horizontal surface are presented. 

Photosynthetic active radiation (PAR) is that part of the solar radiation with a 
wavelength of 0,38 to 0,71 wm, which is used during photosynthesis. The nearest 
meteorological station for which there is data is used and data are presented in Table 2. Due to 
low variation because of low distances in latitude and longitude the territory of Bulgaria does 
not experience big changes in the solar radiation (Climate reference book for Bulgaria, 1978). 
The total sum of PAR is sufficient for the vegetation cycle as it is the most important factor 
for growing lavender that has a direct effect on the quality of the yield. 


Table 2. Hourly and daily sum of total photosynthetic active radiation (PAR) in cloudless sky 
(cal./cm”). Sofia meteorological station 































































































Sunrise | 5.3°- 62° 72: 8.30. 930. 10.2% 113° 12.3. 13.3% 143°. 153° 162° 173% 18.30 Daily 

Month 5.30 6° 7 8° 9? 10.°° as 2° 13° 14° 15° 16” in” 13° 19° sum 
I 2.6 6,4 10,1 13,4 14,8 14,0 11,5 73 2,0 82,1 
I 0,6 5,7 10,8 16,6 20,6 23,0 22,1 18,6 13,3 7,8 1,5 140,6 
1 0,3 4,6 10,5 17,0 22,1 26,2 27,5 27,0 24,8 20,1 12,8 5,9 0,6 199,4 
IV 0,2 2,8 11,0 18,0 23,7 28,5 31,2 32,1 31,6 28,8 23,0 17,6 9,2 1,8 259,5 
Vv 1,2 6,1 13,0 20,0 26,5 31,3 35,0 36,6 35,9 32,2 26,5 18,8 12,7 6,1 0,7 302.6 
VI 2.0 7,5 14,4 21,3 27,5 31,7 34,7 37,2 37,7 33,5 22,4 21,5 15,3 8,0 1,8 316,5 
Vil 1,2 7,6 14,6 20,4 26,4 31,0 35,2 37,0 36,0 32,7 28,1 21,5 14,3 73 1,4 314.7 
Vill 0,2 24 5,9 11,9 18,2 23,5 28,9 32,3 33,5 32,5 29,1 23,9 17,1 9,1 0,4 268,9 
Ix 1,2 4,5 9,1 15,4 21,1 25,8 28,6 29,1 27,9 23,8 17,7 10,0 0,6 214.8 
x 0,2 2,1 5,6 OK) 16,4 20,4 23,2 23,7 22,0 17,1 10,5 4.6 156,7 
XI 0,6 2,7 9,4 14,3 16,8 17,8 15,9 11,9 73 1,6 98,3 
xl 1,4 6,5 10,5 13,2 14,5 12,9 9,5 4.8 0,6 239 

Soil characteristics 


On the territory of the village of Borika there are 11 soil types classified according to 
WRBSR 2014 update 2015 (IUSS Working Group WRB, 2015), with soil texture 
predominating from low to medium sandy clay according to Bulgarian soil texture 
classification (Ne 1,2,3,4,6), followed by Ne 5,7,9 and low sandy-clayey Ne 8 and 10, as well 
as sandy Ne11. The humus horizon's depth varies depending on the soil type (10-30 cm) and 
all the soils are with low depth. According to the depth of the profile the soils with low depth 
(Ne 7,8,9,10,11) and moderate depth (Ne4,5,6) predominate, followed by deep (Ne2,3) and 
highly deep (Nel). The soils in the region have a low texture differentiation (Texture 
coefficient) ranging from 0.9 to 1.5, with the exception of Ne3 and Ne5, which are highly 
differentiated. The soil response ranges from moderately acidic to slightly acidic (pH in HO 
4,8- 6,0). According to the humus content the soils are low to medium humus (1-3,1%), the 
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groundwater level is deep, over 600 cm. Affected by erosion are almost all soil types with the 
exception of (Nel, Ne2, No3). 


Table 3. Matrix of soil characteristics 









































Soil type <0,01 in <0,01 in area under Fens ool Texture | pH Hu ee ale 
Ne} WRBSR 2014, plowed plowed horizon eee a coefficie in mus ae 
u 2015. horizon (%) (%) (om J a ) nt H,0 | (%) 
! | Fluvisols 35 38 30 160 1.3 6 2.1 A 
2 | Phaeozems 34 39 30 150 1.3 5.3 1.2 A 
3 | Planosols 37 60 25 100 2:2 5.5 1.4 A 
4 | Luvisols 39 49 25 95 1.5 5.3 2.4 A 
5 | Planosols 14 51 10 70 4.7 4.8 1.2 B 
6 | Cambisols 40 26 20 65 1.1 5.6 1.3 A 
7 Leptosols 14 14 15 20 1 5.3 1 B 
8 | Cambisols 28 20 20 40 1.3 5.4 1.4 B 
9 Leptosols 16 12 20 40 0.9 5.1 1.7 B 
j Leptosols 25 25 15 15 1 5 3.1 B 
Leptosols 8 8 10 20 1 5.3 1 B 



































Farmers need efficient practical recommendation for decision making and risk 
avoiding. A variant of the limitation method is proposed where a matrix of characteristics is 
created. This is a mountainous region with lands facing natural constraints and there are 
limitations for agriculture. With —A” are represented the soils that can be used for growing 
lavender with a lower risk, corresponding to the FAO S order(S1, $2, S3) and with -B” are 
represented the soils that have limiting factors of significance that make them non-suitable for 
growing lavender, corresponding to FAO N order (N1, N2) for Borika village (FAO, 1976, 




































































2007). 
Table 4. Matrix of land characteristics evaluation 
<0,01 in plowed | <0,01 inareaunder | Humus horizon | Soil profile ese t pH in Humus 
No horizon (%) plowed horizon (%) depth (cm) depth (cm) <0.9 H20 (%) Count Y 
>30 >30 >15 >50 ? >5.2 >I 
ie 
1 Y Y Y y Y ea 7 Tn 
2 Y Y Y x Y Y ¥ Tn 
3 ay x Y 7 N x 7 6t 
4 x Y ¥ x N x Y 6t 
5 N Y N Y N N Y 3ca,ha,t,p 
6 x N Y 7 Y x Y 6cb 
7 N N N N Y Y N 2ca,cb,ha,d,hb 
8 N N Y N Y Y Y 4ca,cb,d 
9 N N Y N N N Y 2ca,cb,d,t,p 
10 N N N N Y N Y 2ca,cb,ha,d,p 
11 N N N N Y Y N 2ca,cb,ha,d,hb 
The small letter represents the soil limitation factors according to the legend: 
n—no limitations; 
ca — clay content in plowed horizon limitations (%); 
cb - clay content under plowed horizon (%); 
ha — humus horizon depth limitations; 
d—soil profile depth limitations; 
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t — texture coefficient limitations; 

p — pH limitations; 

hb — humus content limitations. 

With -¥” are represented the soils that fall under criteria for growing lavender and 
with —N that do not. 

As it is seen, the best soils in the region for growing lavender are 1, 2, 3, 4, and 6. As 
lavender tuft is not growing very deep it can tolerate slight differences in texture coefficient. 
The clay content in number 6 is closer to the limit and it has points close to the maximum sum 
of 7 (26). They all have equal to or more than 5 -¥”’’s. Number 5 is very acidic, with very 
shallow humus horizon and unfavorable soil texture, number 7, 8, 9, 10 and 11 are with very 
low depth and/or unfavorable soil texture and/or humus content and/or acidity. 

For greater visibility in Fig. 1 and Fig. 2 maps of the studied area are also presented. 
The large-scale (1:10000) map of the land (digital archive of ISSAP "N. Pushkarov", digitized 
by EASR, MAF) includes 11 soil mapping units and is shown in Fig. 1. Note that these are 


images not actual maps with scale. Not all of the land is used for agriculture. Some of the 
land is covered with water, forests, roads, buildings. 





Soil type differences Area (ha) 


[ Fluvisols, average depth 51,579 
= 










Phaeozems, average depth 152,433 
Planosols, no erosion 150,693 
Luvisols, low erosion 180,056 
Planosols, low and medium erosion 124,622 
=] Cambisols, medium erosion 11,801 
[__] Leptosols, medium and high erosion 364,628 
[9 Cambisols, medium leached, low and medium erosion 456,568 
an Leptosols, tallow, few and median eralon hens 
eptosols, shallow, medium and high erosion 
1 0 1 2 3 Kilometers 1 Leptosols, low leached, shallow, medium erosion 21,220 
| Ravines and gullies 91,827 


Figure 1. Soil map of Borika village, Ihtiman district, Scale 1:10000 
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Legend 


HB Water area 
HBB Ravines and gullies 
Other territory used by agriculture 
HB Forests 
HB Built-up area 
[_] Meadow 
HB Local road 
[___] Agriculture field 
Pasture ground and shrubs 
1 0 1 2 3 Kilometers SS Agriculture fields' road 
a 








Figure 2. Map of distribution of the land of Borika village according to methods of long-term use 
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Figure 3. Area of soil units by number 
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Conclusion 

It is more difficult to grow lavender in areas facing natural constraints. One of the 
units in the legend are not soils. These are ravines and gullies according to the Bulgarian soil 
map.Most of the area is covered with soil number 10, 8 and 7 (Figure 3). 

The best soils for growing lavender are with number 1,2,3,4, 6 (representing Fluvisols, 
Phaeozems, Planosols, Luvisols and Cambisols) and they cover only 21,1% of the land. They 
are located near the rivers and in specific valleys with small slopes and are most suitable for 
growing lavender. In conditions that face natural constraints growing lavender should be 
cautious. 
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Abstract 

A containers experiment was carried out in four variants with fertilization of ryegrass: 
1th and 2nd trials, with the addition of 1ml and 5ml of organic fertilizer “Siapton” and 3th and 
4th trials: a combination of Siapton (0.5 ml fertilizer in 50 ml of water ) with the mineral 
fertilizer "Kristalon" in quantity of 0.5 and 2 grams respectively. There was used also a 
control. It is accounted a higher accumulation of organic matter (content of C% and N%) in 
leaves of ryegrass after fertilizing in comparison with the control. Accumulation of dry matter 
also is higher in fertilized variants (27.09 + 37.62) compared to the control (25.99g). This 
difference is not confirmed and does not recorded in relation to the roots of ryegrass. 
Acidification of soils is changed towards a higher pH (5.5 + 6.4) in fertilized variants in 
comparison with the control (pH = 5.2). The humus composition was studied by the method 
of Kononova - Belchikova, but apparent differences were not accounted in relation to the use 
of organic and mineral fertilizers. The content of humic acids extracted with pyrophosphate 
calculated according the weight of soil samples is lowest in control (1.02%), where the 
content of fulvic acids is high (0.40%). 


Key words: organic mass, dry weight, humic and fulvic acids 


Introduction 

Ryegrass (Lolium Perennel L.) is a short-lived grassy species. It is used as a forerunner 
of many crops as it develops a powerful root system and thus improves the soil structure, 
increases overall soil moisture and nitrifying ability. Ryegrass has a special place in the cities, 
the landscaping of recreational areas and sports grounds. It also participates in improving the 
aesthetic features of gardens, parks, family houses and amusement facilities. As a one-year- 
old grassy species, ryegrass also has a place into the composition of many fodder mixtures. 
Rich in protein and energy quality high-quality fodder is obtained by cross-sowing of ryegrass 
and white clover (Mihovski and Sabeva, 2011). As a forage culture, it can also be used in a 
number of agroforestry systems (Kachova et al., 2018). And this: increasing the yield and 
quality of forage and creating weed-free areas is an important part of the tasks in modern 
agriculture. Ryegrass is also used for afforestation of lead-contaminated lands to cover the 
objectives of phytoremediation (Bgantzova, 2009). Petrova (1992) uses the fertilization of 
ryegrass during its cultivation on heavy metals contaminated substrates of phosphogypsum 
and pyrite sludge in order to improve its state. Ryegras reacted violently after the fertilizers 
were brought into. Fertilization has a positive effect on its state. In general, it can be affirmed 
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that the use of fertilizers in ryegrass cultivation leads to an increase in yield, improvement of 
its condition and reduction of negative effects (Vassilev, 2004, 2006, Dimitrova, 1995). 
Fertilization can also improve soil properties (Valchovski et al., 2014). 

The aim of the present study is to trace the effect of fertilization of ryegrass with 
organic and mineral fertilizers on the accumulation of biomass and the soil properties. 

For this purpose organic fertilizer are used - "Sapton", which is new for the Bulgarian 
market and mineral fertilizer "Kristalon", rich in N, P, K and other nutrients. 


Material and Methods 

In the spring of 2017, ryegrass (Lolium Perennel L.) was seeded in containers of 
500ml volume in a dose of 0.2 grams seeds per container. The soil was sieved and cleaned up 
by roots. Soil type is Vertisols and it was taken from the yard of the Forest Research Institute 
- BAS. After 29 days fertilizers were applied and after 30 days treatment was repeated. The 
following experimental variants were used: 

- Treatment]. Fertilizing with 1 ml of "Siapton" organic fertilizer dissolved in 50 ml of 
water. 

- Treatment 2. Fertilizing with 5 ml of "Siapton" organic fertilizer dissolved in 50 ml 
of water. 

- Treatment 3. Fertilizing with 0.5ml organic fertilizer "Siapton", dissolved in 50ml of 
water and 0.5g ,,Kristalone” mineral fertilizer. 

- Treatment. Fertilizing with 0.5ml organic fertilizer "Siapton", dissolved in 50ml of 
water and 2g “Kristalone” mineral fertilizer. 

- Treatment 5 - control without fertilizing. 

Each treatment was repeated in 7 repetitions. The control was in 14 repetitions. Thus, 
all containers with the control were 49 in numbers. 

Throughout the experiment, constant moisture was maintained, with even watering of 
all containers. The organic fertilizer "Siapton" was developed on the basis of naturally 
hydrated proteins from firm "Chemtura Europe Limited", UK. This organic fertilizer is 
specialized in increasing yield and quality of production and activating protection forces of 
crops. The composition given by the manufacturer is: total organic ammonium nitrogen up to 
4%; organic carbon - 25%; amino acids with animal origin - 54.34%. The “Kristalon” mineral 
fertilizer is a production of firm "Hydro Agro Rotterdam", the Netherlands. Contains N, P, K, 
Mg as the main ingredient and other trace elements. 

In September, the ryegrass was removed with the roots, washed and dried at 105 ° C to 
a constant dry weight. Plant samples (roots and leaves) for C and N content were taken. 
Carbon was determined with a bichromic mixture (0.4N K2Gr207 + k.H2SO4) in a plant 
sample of 0.05g through the method for determination of carbon in soil samples (Donov et al., 
1974). The nitrogen in the plant samples was burnt with k.H2SO,4 and H2O2 and determined 
with a Parnas-Wagner apparatus with 50% NaOH and titrated with 0.02N H2SO4. Soil 
samples were analyzed for C - Thurin method, N - Keldahl method, and composition of 
humus - by the method of Kononova - Belchikova (1961). 
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Results and Discussion 

After dry weight measurement, it has been found that fertilization produces a 
corresponding effect in the accumulation of biomass in the grassland. The results are 
presented in Fig. 1. On the figure are also given the weight of the roots. 


Dry weigh of roots and leaves of ryegrass 
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Figure 1. Dry weigh of roots and leaves of ryegrass 


The biomass of the roots far exceeds that which is accumulated in the leaves. It 
confirms the fact that ryegrass develops a powerful root system. The lowest is the weight in 
the unfertiliezed variant - the control. The best is accumulated leaf biomass in the ryegrass of 
Treatment 2, in which we used the largest amount of organic fertilizer "Sapton" - 5 ml, 
dissolved in 50 ml of water. In the roots, this effect is not observed. Fertilization, including 
with mineral fertilizer, does not affect the accumulation of biomass in the roots. On the 
contrary, the untreated variant (the control) also accumulates a large amount of root mass. 

Particularly well fertilization influences the content of carbon and nitrogen in the 

leaves - fig. 2, especially that of carbon. 
The control has the lowest amount of carbon stored in the leaf mass of ryegrass. Fertilizing 
with “Siapton” organic fertilizer improves the carbon and nitrogen accumulation in the leaves, 
especially when using more fertilizer - Treatment 2 (Sml). Especially advantageous is the use 
of “Kristalon“ mineral fertilizer (Treatment 3 and Treatment 4). The maximum amount of 
carbon and nitrogen we have at the experiment with the lower dose of “Kristalon“ (Treatment 
3). Perhaps the higher dose is not well combined with the use of the organic fertilizer 
"Siapton". 
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This is not the case with the accumulation of carbon and nitrogen in the roots of 
ryegrass - fig. 3. Even on the contrary, the control also has significant carbon and nitrogen 
content in the roots. Here we found a lack of effect from fertilization with here used specific 
fertilizers with specific doses on biomass accumulation and carbon and nitrogen content in the 
ryegrass roots. 
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Figure 2. Percentage of carbon and nitrogen in the leaves of ryegrass 


In soil studies the main results are given in Table 1. The results are mean values of the 
7 repetitions. 
Table 1. Main characteristics of soil 




















Experiment pH C% N% C/N 
1 135 4.6 0.5 9.2 
Z 6.14 3.9 0.4 9.8 
3 6.34 3.9 0.4 9.8 
4 6.29 338 0.4 9.5 
K 5.24 3.8 0.4 9.5 




















Soil acidity is favorably influenced by fertilization with organic and mineral fertilizers. 
From a moderate acidic reaction in the control (pH = 5.24), the acidity decreases to a slightly 
acidic (Treatments 2, 3, 4) which acidity is more favorable for vegetation development. The 
carbon and nitrogen content of the soil is also affected by fertilization. Particularly beneficial 
is their content in Treatment 1. But the t-test done did not show statistical proven differences 
between contents of C and N in fertilized variants in relation to the control. Overall, however, 
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in all Treatments, the C/N ratio is very low (<16), indicating that the mineralization processes 
in this soil are delayed and the humus formation proceeds at a low speed. 

The analysis of humus composition according Cononova - Balchikova was also made. 
It determines the total content of humic and fulvic acids using a mixed solution of 0.1 M 
NasP207 and 0.1M NaOH (Table 2). 
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Figure 3. Percentage of carbon and nitrogen in the roots of ryegrass 
The table shows the percentages of the total carbon extracted with the pyrophosphate 
extract; the amount of humic and fulvic acids; “aggressive” fulvic acids and carbon in the 


residue. Each percentage is calculated towards the weight of soil sample, taken for analysis. 


Table 2. Humus composition 




















Experiment Cin Humic acids Fulvic Cinthe | “Aggressive” 
pyrophosphate % acids residue Fulvic acids 
extract % % % 
% 

1 1.59 1,32 0.27 2.99 0.20 

2 1.72 1.25 0.47 1.66 0.19 

3 1.62 1.08 0.54 2.27 0.08 

4 1.42 LAS O27 2.40 0.23 

K 1.42 1.02 0.4 2.4 0.08 
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Carbon in the pyrophosphate extract is low (~ 1%). On the contrary, the non-extracted 
carbon retained in the residue is in a higher percentage (~ 2%). Fertilization favorably affects 
the content of humic acids. Humic acids are more stable organic acids in composition to 
organic matter. They sequester of larger amounts of carbon in its composition. In this sense 
the effect of fertilisation is in relation with sequestration of carbon in the system.The lowest 
is the percentage of humic acids in the unfertilized variant (control). On the contrary, the 
percentage of the more unstable and more easily mobile fulvic acids in the soil is high in the 
control (0.4%). Aggressive fulvic acids, however, are not high in this non-fertilized variant, 
probably due to the low level of carbon extraction in the pyrophosphate extract. 


Content of humic acids in soils calculated as 
percentage from total carbon in soils 
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Figure 4. Content of humic acids in soils calculated as percentage from total carbon 
in soils 


If, however, we calculate the amount of humic acids not in relation to the quantity of 
the sample but as a percentage of the total carbon content in the soil (Table 1), a clear picture 
is obtained - Fig. 4. 

Clearly is hilighted the result, that the amount of humic acids is the lowest in the non- 
fertilized variant - the control. Hence, we can conclude that fertilization stabilizes organic 
matter in the soil, improving the amount of more stable forms of carbon - the humic acids. 
Sequestration of carbon is increased. 
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Conclusion 

The use of organic fertilizer "Siapton" alone and in combination with "Kristalon" 
mineral fertilizer has a beneficial effect on the accumulation of biomass in the ryegrass, as 
well as on the properties of the soil. Fertilization results in: 

- the accumulation of more dry matter in the leaves; 

- higher carbon and nitrogen content in the foliage; 

- improving soil acidity; 

- better store of organic matter; 

- a more favorable humus composition in the soil with a higher content of the humic 
acid fraction. 

No repercussions of fertilization on the accumulation of biomass in the roots of 
ryegrass are detected. No statistical proofs were found. 
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A meeting related to the European scientific program COST (European Co-operation 
in Science & Technology) took place from 25" and 27" September, 2018 r. in the Congress 
Center of Flamingo Grand Hotel (Snapshot 1). Representatives of 39 countries and 
organizations as the European Forest Institute and the United Nations Economic Commission 
for Europe/ FAO participate in the cost action PESFOR-W, Snapshot 2. Chairpersons of the 
action are D-r Gregory Valatin (United Kingdom), and Vice Chairmen — Prof. Gebhard 
Shueler (Germany) [1]. The hosts of the meeting were Prof. Dr. Mariana Lyubenova from the 
Faculty of Biology, Sofia University and Assoc. Prof. Alexander Shikalanov from the Faculty 
of Information Sciences, University of Library Studies and Information Technologies, who 
are the members of Management Committee of the Action. 





* Sw, 
eri Bas = 
Snapshot 1. The meeting in action. Snapshot 2. Participants infront of the Green Center. 


The main objective of international scientific cooperation is to collect and validate 
good practices of forest ecosystem services assessment for the quantity and the quality of 
water in adjacent areas (selection of complex water flow and water quality indicators). The 
expected practical outcome of international science cooperation will be to synthesize 
knowledge, promote joint research to improve Europe's capacity to use Payments for 
Ecosystem Services (PES) for achieving the Water Frame Directive (WFD) objectives and 
other policy goals through incentives for planting forests to reduce diffuse pollution from 
agriculture. Also the result will be the development of PES Schemes, the creation of a 
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protocol and algorithm for calculating the cost- effectiveness of forest systems and the 
schemes developed. 

The representatives of Bulgarian institutions presented their original achievements 
related to the management of water resources (Basin Directorate "Black Sea Region"), 
development of payment schemes (Institute for Economic Research, BAS and WWF- 
Bulgaria), management of forest resources in Bulgaria, mapping and assessment of forest 
ecosystem services and legislation developments (Executive Forest Agency, Forest Institute, 
BAS). The first results from the Bulgarian project ,,Development of an ecosystem services 
assessment scheme, their effectiveness for purification and protection of water and other 
natural components in the regions“ [2] (National |Scientific Fund, Ministry of Education and 
Sciences), were presented (Sofia University, University of Library Sciences and Tecknology). 
Some of the first dendromers in natural forests were tuned, installed and put into action, that 
monitor the 24-hour radial growth of Quercus frainetto Ten. and Quercus cerris L. Long-term 
observations would make it possible to create a database and participate in international 
projects (Snapshot 3 - 4). Professor Margaret Shannon of Baldy SUNY Buffalo Law School 
presented the subject "Critical Transformational Deliberative Science: A Critical Element for 
PESFOR - F for W“. 


~t i 
a Sees 





Snapshot 3. Dendrometer on Q. frainetto tree. Snapshot 4. Dendrometer on Q. cerris tree. 


One of the focus of meeting in Albena is the sistem of field protective forest belts in 
Dobrudja [3-5]. The role of wetlands for water regime of the regions, for the biodiversity 
protection, as a bio-carbon storage and for the natural components purification, is also on the 
focus of international scientific cooperation. The Bulgarian experience in the construction, 
maintenance and management of field protective forest belts in the Bulgarian part of Dobrudja 
was reported (Northeastern State Enterprise, Forest Protection Station -Varna, State Hunting 
Farm —Balchik, State Forestry Farm - General Toshevo). The system of belts in Dobrudja as a 
unique phenomenon in the agroforestry practice and an original ecology-biological 
phenomenon caused a great interest of participants (Snapshot 5). Everyone agreed that the 
forest belt system in Dobrudja is not only a national wealth, because of its uniqueness and the 
complex of benefits (services) it offers. This national wealth must be preserved, maintained 
and developed. It is necessary to develop PES schemes that provide funding for their 
sustainable management, as well as for the reconstruction and creation of new belts. The high 
professionalism and experience of Bulgarian foresters needs to be exported to Europe as a 
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good and successful practice that will greatly help to create protection belts in the other 
European countries and to find funds supporting the maintenance of Bulgarian shelter belts. 

The biodiversity of wetlands in the Shabla and Balchik Municipalities as one of the 
main factors determine ecosystem sustainability in the context of climate change as well as 
the ecosystem services provided by wetlands were also disscussed on the field trip. 





Snapshot 5. Disscussion on structure and lenefits | Snapshot 6. Duranculak Lake and its biodiversity. 


of forest belts. 
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Snapshot 7. Longozes in Baltata Manage Reserve. Snapshot 8. Archeological findings on B. Island 


The protected area "Durankulak lake" (Snapshot 8), which is one of the most significant and 
best preserved coastal wetland in Bulgaria, registered under the Ramsar Convention and 
included in the Natura 2000 network was visited [6]. It is a natural liman fed by karst waters 
and separated from the Black sea by sand dunes and a beach strip. The Via Pontica road 
passes through and making it a suitable destination for ornithological turizm. Participants in 
the international meeting visited also the managed reserve "Baltata", located in close 
proximity to the Albena resort, where the northernmost longoze forests occupy a place 
(Snapshot 7). These forests are unique communities in the vegetation of Bulgaria with 
significant biodiversity and richness of conservation - significant species [7]. The forests are 
formed in the lower course of Batova River, which originates from the village of Kumanovo, 
the Frangen Plateau. The total length of Batova River is 39 km, the minimum run-off is 200- 
300 1/s. Longozes occurred at the end of the glacial period on quaternary deposits, rich, wet 
and deep soils. An important factor for their normal development is the soil moisture. The 
presence of evergreen and warm-loving species, as well as a more diverse species 
composition of tree layer, also a greater amount of wood and grass lianas distinguish these 
forests from the other riparian forests in Bulgaria (Regional Inspectorate of Envinroment and 
Water — Varna, Shabla Municipality, Green Educational Ccenter). The managed reserve was 
commented on, as a natural example for the role of forests for the surface water purification. 


Bulgarian Journal of Soil Science® 2018 Volume 3. Issue 2. www.bsss.bg 152 


International Meeting 





There are 5 villages in the course of the river. The town of Balchik and v. Kranevo are resorts 
and Albena is a preferred resort on the seashore. They are located in the lower reaches of the 
river. During the summer season, the recreational load is strongly expressed. The agricultural 
lands also occupy a large area. Two roads from the State road network pass through the valley 
of the river. Extremely high pollution of the surface water flowing into the Batova river bed 
can be expected. According to the data from the national monitoring, the Batova River valley 
is characterized by medium to slightly polluted water and pure water flows into the sea. In this 
case, for the good ecological status of surface water in the Batova basin and coastal seawater, 
the regulating and supporting ecosystem services of longoze forests are of great importance. 
The maintained reserve "Baltata" also contributes significantly to the unique natural 
conditions that are offered to the visitors - a unique combination of preserved nature and a 
modern tourist resort providing various forms of recreation and ecotourism. 

During the field trip, the rich cultural and historical heritage of Dobrudja and 
Durankulak as a complex archaeological site, where the highlights of Dobrudja's prehistory 
have been presented - an attractive center for the cultural and archaeological tourism [8-10]. 
The Big Island in Durankulak Lake and the Museum Collection in the Green Educational 
Center — Shabla were visited. On the Eneolithic settlement mound (Big Island), the oldest 
stone architecture in Europe (Dr. I. Weissov, NAM, BAS) was found (Snapshot 8). The 
activities of Green Educational Center at Shabla Municipality for the environmental education 
is the most direct route for PES schemes realization and for sustainable development of 
region. 
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